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Summary—Experimental data are given for the transfer of material by diffusion from a continuous 
phase to dispersed drops. The data are given for three conditions: (a) during formation of the drop, 
(b) for transfer to the drops rising through still liquor, and (c) for the transfer to the drops rising counter 
to the downflow of the continuous phase. The transfer effected during formation is almost independent 
of the time of formation of the drop. The overall transfer coefficient K (based on the average area of 
interface during formation) decreases with increasing time of formation of the drop but is practically 
independent of the size of the drop. Smaller drops attain more nearly to equilibrium with the liquor 
because of the increased area per unit volume. 

The value of A during the rise of the drops increases with drop size and is directly influenced by the 
relative velocity of the drops and the continuous phase. Comparison of the results with those anticipated, 
assuming that the drops rise as spheres or that they renew their interface each time they rise through 
one diameter, suggests that the real conditions lie between these extremes. The time of exposure of the 
interface is a main feature in determining the actual value of K. The results are applied to the case of 
a perforated plate extraction column. 


Résumé@— Données expérimentales sur l'échange materiel par diffusion entre une phase liquide continue 
et des gouttes de liquide qui y sont dispersées. Ces résultats correspondent A trois cas distincts: 

a) Formation de la goutte, 

b) ascension de la goutte dans le liquide immobile, 

¢) ascension de la goutte & contre-courant de la phase continue descendante. 
a)—Durant la formation de la goutte, léchange est A peu prés indépendant de la durée de formation de 
la goutte. Le coefficient global de transfert A (basé sur la valeur moyenne de la surface de contact 
durant la formation) diminue avec une durée croissante de formation mais reste pratiquement indépen- 
dant de la taille de la goutte. L’équilibre avec la phase continue est atteint plus rapidement avec de 
petites gouttes A cause de l'accroissement de surface par unité de volume. 
b)—La valeur de K croit avec la taille de la goutte. 
c)—La valeur de K est en proportion directe de la vitesse relative des gouttes et de la phase continue. 

Les résultats expérimentaux paraissent indiquer que l’échange intervient dans des conditions inter- 
médiaires entre les deux extrémes: ascension de la goutte conservant sa forme sphérique ou renouvelle- 
ment complet de la surface de contact chaque fois que la hauteur d’ascension atteint un diamétre. C'est 
la durée de contact & la surface de séparation qui est le facteur déterminant pour la valeur actuelle de 
K. Les résultats sont appliqués & une colonne d’extraction a plateaux perforés. 


INTRODUCTION AND Scope oF PAPER Eq. (2) avoids the determination of the concentra- 
There have been a number of papers describing the tion at the interface by the use of overall transfer 


performance of small laboratory liquid-liquid ex- coefficients A, which are related to the film coef- 


traction columns [1], [2], [3], [6], [8], [9], in which ficients by the equation, 

the general treatment has been to apply the two ' ' H : 
film theory of mass transfer by diffusion, and to ( 9 , he « %, 
express the transfer obtained in the form of a general 


(3) 
rate equation such as: For any apparatus such as a column eq. (2) is 
Nit = k(Cy— Cad 8S = by (Chi — Oy). (1) 

On the assumption that equilibrium conditions 
exist at the interface and where the equilibrium where AC,, is the log mean driving force at the 
relation is a straight line this relation may be written, jnlet and exit of the unit. Since S is often very 
Nit = Ky(Cq — CI) S = Ky (CF — Cy) 8. (2) difficult to measure it is frequently replaced by the 

* Present address: Monsanto Chemicals, England. term a.V, where a is the interfacial surface per 


conveniently written in the form 


Nt=K.S.AC,, (4) 
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unit volume of the column and V is the volume of 
the unit. This gives the general transfer equation as, 


Nt=K.a.V. ; (5) 


This equation is widely used to calculate values 
of Ka from experimental data but as both A and a 
difficult to 
the results are not easy to interpret. 


are variables and a is very measure, 

The disadvantage of the above equation is that 
it assumes steady state conditions whereas in practice 
one liquid is usually dispersed within the other and 
the concentrations and rates of mass transfer may 
vary with time. The simplest case of unsteady state 
mass transfer is the absorption of a pure gas in a 
still liquid and this case was studied by Hieste [4)}. 
When the gas and liquid are first brought into 
contact there is a very high rate of absorption which 
the 


depth of liquid to be infinite and Granams law of 


progressively decreases with time. Assuming 


diffusion to apply 


(6) 


Integrating we obtain 
yt 
; 2(C" ' | Dt, 


7 


amount absorbed in time ¢ and 


ky 2| ¥ (7) 


where / 


nt, 
The above differential equation can also be solved 
for various simple geometric shapes and such solu- 
NewMan |12). The 


equations are important since they show & to be a 


tions have been given by 
function of the time of contact or exposure of the 
interface. They may be applied to extraction if it 
is borne in mind that there are two films instead of 
one but the process is often complicated by turbul- 
ence. 

There is some evidence in the literature to suggest 
that K varies with the time of exposure of the inter- 
face but no quantitative information, and the present 
work was undertaken in order to study the influence 
of this variable and to obtain a more complete 
In the 


first paper the mechanism of mass transfer into 


knowledge of the mechanism of extraction. 


drops is considered while the second is concerned 


with extraction between continuous liquid streams. 


Part 1. Mass transfer into single dispersed drops 
The general equation for counter current extraction 
columns does not take into account the fact that the 
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extraction takes place in three distinct stages. The 
first stage occurs when the drops form, the second 
when they rise or fall through the continuous phase, 
and the third when they coalesce. This cycle may 
occur only once as in spray towers, or be repeated 
as in sieve plate towers. The present work was 
aimed at separating these stages and studying them 
individually, data for the first two stages being 
obtained by direct experiment, and that for the 
third stage by difference. 

Very few papers have been published on the 
transfer to single drops and the results are not as 
yet conclusive. SHerRwoop [7], passed drops of 
benzene or methyl isobutyl ketone containing a 
small percentage of acetic acid as solute into a column 
of water and measured the change in concentration 
of the inlet and outlet streams. He expressed the 
overall transfer coefficient in the form of eq. (4), 
using the change in concentration multiplied by the 
flow rate for NV, and the total surface of the drops 
in the column at any one time for 8S. He found that 
A increased with drop size, and with the time of 
formation of the drop. By plotting the ratio of the 
unextracted solute to the total solute it is possible 
to extract, against the column height he found that 
some 40% of the extraction measured occurred during 
the formation of the drops. West [ll], using a 
similar system and apparatus, however, obtained 
results which differed very considerably from those 
of Suerwoop. Licut and Conway [5] allowed drops 
of water containing acetic acid to fall through a 
small column of isopropyl ether, methyl isobutyl 
ketone, or ethyl acetate. Using different heights of 
column they extrapolated back their data to show 
that 18-20% of the total extraction occurred during 


drop formation. 

The present work has been arranged to study 
these features individually and is conveniently ar- 
ranged in three series: 


Series 1. Mass transfer into forming drops. 
2. Mass transfer into drops rising in a still 
continuous phase. 
3. Mass transfer into drops rising in a con- 
tinuous phase moving against the drops. 


Two systems have been used, benzene-benzoic 
acid - water, and benzene -propionic acid - water, the 
76 respec- 
tively for the concentrations used in the experiments. 


Throughout benzene was used as the dispersed phase. 


distribution coefficients being 0-068 and 2:- 
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A few experiments have been made with methyl 
isobutyl ketone-benzoic acid-water where the distribu- 
tion coefficient is 0-014. 

The results obtained have been used to examine 
some previously reported data on perforated plate 
and spray towers. 


EXPERIMENTAL WorK 
Series 1. Transfer to forming drops 

For precise measurement of the extraction occuring 
during drop formation, it would be desirable to 
analyse the drop immediately after it is formed. This 
is very difficult and the method adopted has been 
to form a drop on a jet in the continuous phase, and 
as soon as it is formed to push it out of the system 
Two forms of equipment have been used as shown 
in Fig. | and 2, the first for conditions of slow forma- 
tion of drops and the second for more rapid formation. 

As shown in Fig. 2, the benzene phase flows 
from a small reservoir through a needle valve to 
form a drop at the end of the glass jet fixed in the 
hase of the column holding the continuous phase. 
A side tube of 4 mm i.d. is sealed just below the jet, 
so that if the two way cock is turned the flow of 
benzene is cut off and the drop is forced out of the 
system by the hydrostatic head of the continuous 
phase. The procedure was to set the needle valve 
to give the desired rate of drop formation and when 
the drop had reached the desired size as shown by 
a scale at the back of the tube, the cock was turned 
and the drop pushed out of the apparatus and collected 
in the burette as shown. This cycle was repeated 
until sufficient liquor had been collected for analysis, 
the liquor in the column being stirred after every 
ten drops. The drop volume was obtained from the 
number of drops and the total volume and time of 
This method was 
only suitable for times of formation greater than 


drop formation were recorded. 


3 sec. 

For more rapid drop formation the automatic 
apparatus (Fig. 2), was developed. Here the head 
of liquor in the column is balanced by an equivalent 
head on the inlet and outlet tubes. If the pressure 
at the top of the column is slightly reduced a drop 
will form on the jet, and when the pressure is increased 
again the drop will be forced out of the column. By 
the use of two ball valves the drop is drawn up the 
inlet tube, but is pushed back out of the side tube 
to the measuring burette. Drops were formed and 
removed by fluctuating the pressure at the top of 
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the column using a small piston pump as shown. 
The pump was driven by a geared motor so that a 





Apparatus for study of extraction into forming 


Fig. 1. 
drops. 


Na lech 
| Gor eching 
burette 





————~ 
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fo 
Storege bottle 





Fig. 2. Apparatus for study of rapidly forming drops. 


wide range of rates of drop formation could be 
obtained (j to 3 secs). By fitting a variable throw 
eccentric the actual drop volume could also be 


varied. 
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Four sets of experiments Series | A-D_ were 


made at constant drop size, and varying rates of 


drop formation, and one set 1 E, with constant drop 


rate but varying drop size. In these tests benzoic 
acid was extracted from the water phase to the 
dispersed benzene drops. A further set 1 F was made 
using propionic acid as the solute. The analysis in 
all cases was accomplished by shaking the benzene 
sample with excess NaOH, and titrating the mixture 
with HCl, using thymol blue indicator. In all cases 
a blank with the 


quantity of HCl and NaOH, and the benzene nor- 


titration was carried out same 
mality was obtained from the difference between 


the two titres. 


Methods of calculation 
The two quantities of greatest interest are 

1. The total extraction occuring during drop 
formation and withdrawal. 


2. The rate at which the extraction takes place 


extraction—This is conveniently 


1. The 
expressed as a percentage of the extraction that 


amount of 


would be obtained if the drop attained equilibrium 


with the surrounding liquid. 


o. Equilibrium Final normality « LOO 

Equilibrium normality 
where the normality of the collected drops is used 
as a measure of the concentration. 

In the operation of the ap- 
paratus there is a small volume 
r above the withdrawal point 
which is filled with the dispersed 


liquid during drop formation. 





This liquid, however, does not 
enter the drop but is pushed out 
when the drop is collapsed. If V be the volume of the 
actual drop formed, then in each cycle a volume 
(Vv) of liquor enters the apparatus from the 
dispersed liquor reservoir. This is the volume of liquor 
that is actually pushed out and measured in the 
burette, so that the percentage approach to equili- 
brium based on the volume Vo of the actual drop 
is FE, 
Final normality (V+ rv). 100 (9) 
Equilibrium normality 1 
It is convenient to express all results for a constant 
concentration of the aqueous phase. H/, for the system 
varies with concentra- 


benzene-benzoic acid-water 
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tion, and although for any one experiment the 
concentration change is very small, it is necessary 
to allow for the change in H throughout the ex- 
periments in any series as they were not carried out 
at precisely the same concentration. 

For this system H A, is small compared with 


lL ky, 


dent on the physical properties of the system and 


and we may take K, k,. But &, is depen- 


may be taken as constant over the concentration 
range considered. Hence we may take A,, as constant. 
Total x kK, 


phase). For dilute concentrations this driving force 


extraction (driving force in water 
may be taken as equal to the actual concentration 


in the water phase. 
Hence, 
Final benzene normality x Water phase normality (10) 


since the benzene normality represents the mass 
transferred, 

The results in these series have been corrected 
to a water phase concentration of O-O11, for which 
value H equals 0-162. The corrected approach to 
equilibrium for these standard conditions will be: 


Final benzene normality (V+ vr) 100. O-OL] 
0-162 Water normality * 


E (11) 

This correction is very small and hardly necessary 
for forming drops, but is of greater importance for 
the experiments with rising drops. With propionic 
acid as solute H is sensibly constant and no similar 


correction is required, 


2. The 
interfacial area is continuously changing, but an 


rate of extraction—With forming drops the 
average value may be calculated in the following 
way. First we must find an average time of e.- 
posure ¢, for the drop that is formed and collapsed 
in a time ¢,, the time of formation being equal to 
the time of collapse. Assume further that the volume 
of the drop increases uniformly with time and that 


the drop is spherical. 

Let R be the radius of the drop at a time }t,, 
when the drop is fully formed. 

Let M be the volumetric rate of flow, so that the 
volume at time ¢ is M.t. 

Let r be the radius at any time f, so that M.t = 


2 
‘ a. and dt S28 ee 


The interfacial surface at time ¢, S 
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The average surface between 0 and ¢, is given by, 


ty 
"4.2. .di 
a 


dr 


ii 


The actual time of formation and collapse of the drop 
t, is given by 

}t, (13) 
> 

where ¢, is the measured time of one cycle during 


which a volume (Vv) of liquor flows into and 


out of the unit. 


If in any interval of time ¢ to ¢ + dt, a surface 
ds is formed then it will be exposed for a period 
‘, —¢ during the cycle. Hence the average time of 


exposure ¢, is given by 


Noting that VM. that (t,- ¢) 


4.7.(R® — r') 


s 
thds = f(t 


whence ~ t, (1A) 


(Sa.r.di 
0 


ic. the average time of exposure is { the time of 


formation and collapse. For the automatic machine 
the time of formation was equal to the time of collapse 
but this was not accurately so with the hand apparatus. 
Eq. (14) will, however, apply in both cases. 

coefficient 
of the 


average interfacial area exposed per unit time. The 


the 
expressed in 


Using these relations average 


of transfer can be terms 


are evaluated 


data for one run with benzoic acid 


below. 
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Example calculation 


Measured time for formation and collapse of 


2? secs 


drop ¢,. 
Drop volume V+ v 0-085 ml 
Actual drop volume V 0-067 ml 
Final normality of drop sample 0-00053 


Water phase normality . O-OLIG 


O-085 


0-067 


OMS 
0162 


ool 1a 
OOLL6 


o Equilibrium £ 


O39, 


on . O67 
Time for formation and collapse of drop — 2» . 
OUSD 


l “DS SeCs, 


Final radius of drop R, O67 0-021 ems. 


12 2 
*- ; t- = 048 sq cms, 


Average surface exposed S 
Since a normal solution of benzoic acid has 1 gm mole litre. 


the moles of acid transferred to the dispersed phase 


0-085 O00053 


Loo 
3600 
0-48 


214» 10 


0-085 O00053 


Rate of transfer 


1-58 


-om mols hrem?’. 


The water phase concentration does not alter and may 
be taken as the driving force since the benzene phase initially 
had no benzoic acid. 


Driving force O-OLL6 © 10-8 gm mols ce. 


ae . ym mols 
Hence the average transfer coefficient K,, 


0-214 » 
OOLI6 


hrom? 
lo? 


gm mol ce lo * 

18-4 gm mols hr cm? (gm mol ce). 
Giving K,. — 0-61 lb mol hr ft? (lb mol ft*). 
Experimental results for transfer to forming drops 
The experimental results are expressed by consider- 
ing the influence of, time of drop formation and 
collapse t,, the drop size, and the physical properties 
of the system on the rate of transfer and on the 
amount of transfer. A summary of the experimental 
data is given in Table 1. 
The 


transfer effected E is shown plotted against ¢, for 


/. Ej fe ct of time of drop formation and colla pse tl, 
the two systems in Fig.3 and 4. The points for 
propionic acid lie on a smooth curve and although 
those for benzoic acid are more scattered a similar 
It is seen that 
the 


curve has been drawn through them. 
E increases with ¢,, but for values of ¢,>5, 
increase is very small. £ is numerically small for 
both systems, and since the total extraction increases 
very slowly with ¢,, the rate of extraction must 
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Table 1. Benzene-benzoic acid-water. 





Measured 
drop 


rolume 


mi 


O-0685 
0-069 
0071 
0-060 
0-060 
0-068 
(0887 
(0835 
O-OS845 
00-0846 
(-085 
C-OS864 
(1-099 * 
1090 
LIS 
m120 
123 
14 
120 
1S! 
127 
134" 
136" 
138" 
(et 
0-064 
O74 
0-08 
One 
102 
0-123 


mil 


O-0505 
0-051 
005: 
O51 
O-O51 
0-050 
O-O707 
O-0655 
0-665 
00-0666 
0-067 
O-0684 
-OS1 
o-oo! 
0-100 
im 1o2 
105 
0-006 
mil 
mils 
109 
116 
11S 
mize 
0-036 
46 
0-056 
One. 
mos 
O-084 


Average 


exposed 


interface 


. 
cm 


305 
O-4 
41 
O-3o 
O40 
39 
5 
O47 
48 
(48 
(48 
49 
(4 
O59 
62 
0-63 
i 
O-ol 
67 
675 
O06 
O-6S 
One 
7 
O-3l 
O37 
42 
455 
wl 
5D 


Measured 
overall 
drop time 


Actual 
overall 


drop time 


aecn 


Mean 
water 


normality 
<10-2 


1-095 
VW 
“16 
‘ll 
VW 
“16 


Final 
benzene 
normality 

10-4 


Equi- 


librium 


(gm mols cc) 


0-33 
0-37 
0-36 
0-40 
38 
(44 
O33 
O34 
0-36 
O37 
O-39 
415 
625 
(-334 
O3 
326 
37 
39 
(34 
38 
39 
38 
514 
(68 
0-46 
4 
(44 
4 
(376 
(+366 
Ose 


gm mols 
hr om? 


K,, 


lbh mola hr jt 
lbh mols jt® 


1-46 
1-28 
OO] 
OSS 
O67 
O45 
1-38 
1-19 
a2 
735 
(62 
43 
212 
1-2 
OOS 
72 
O63 
Odi 
76 
O-66 
45 
OSS 
192 
140 
a4 
isl 
(SO 
(ESS 
o~79 
Tye 


S25 





Table 1 (continued), Benzene-propionic acid-water 





Measured 
di op 


rol mame 


mel 


(Pray 
0-089 
Oe 
Tete 
0-093 
1-006 * 
0-015 * 
0-096 * 


ietual 
drop 


vol nme 


mil 


OTIS 
mol 
072 
O74 
OTS 
078 
077 
0-078 


Areradge 


‘ Pr poat d 


interlace 


om 


ON 
405 
Os 


Measured 
overall 


drop time 


arcs 


765 
1-43 
1-44 
1-08 
3-03 
424 
6-75 


1O-S 


Ietual 
overall 


di op ferme 


wren 


61s 
O-s82 
1-15 
1-6 
244 
344 
546 
SS 


Mean 
water 


nor mality 


Final 
henzene 
normality 

jo? 


6-73 

74 

| 
10-05 
1S 
12-6 
12-0 


16-3 


Kquilibre ame 


10-65 


qm mols 


hreom 


Kp 


lh mols hr it? 
lh mols 1? 


(gm mols cc) 


40 
$2 


27-6 


131 
145 
OOS 





* Indicates experiments with hand apparatius 
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1952 
decrease very rapidly. This is shown quantitatively 
in Fig. 5 and 6, where the transfer coefficients K,, 
for benzoic acid and Ky for propionic acid are shown 
plotted against t,. The points lie on reasonably good 
straight lines which may be expressed as: 


K, = 0-89 1 5" Benzene-benzoic acid-water 


Kp 


0-96 1/4?7 Benzene-propionic acid-water 














§ 6 ’ é 9 


Time of drop formation and withdrawa 


ene 


Fig. 3. Variation of extraction with time of drop forma- 
tion (for drops of benzene forming in 0-011 N benzoic 
0-07 ml; AO-114 ml. 


acid). Drop Volume 

It is of interest to compare these results with 
the Hiasie theory which as mentioned was developed 
for the absorption of a pure gas into an infinite depth 
of liquid across a plane interface. According to this 


theory 
(15) 
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Using a paper by WiLke [10], and assuming P,/D, 


a constant, for the solvents benzene and water, 


values of D have been calculated for 15° C. 

2-69 « 1075 ft?/hr 
D py 5-28 x 10 5 ft? hr 
D Propionic acid in water 3-4 x 10~° ft®/hr 


wp 
Dp» 6-66 « 107° ft®/hr 


D Benzoic acid in water 


wh 
Benzoic acid in benzene 


Propionic acid in benzene 














. 5 6 4 8 9 
l me of Orog rmotion ond withdrawal 
Fig. 4. Variation of extraction with time of drop forma 
tion (for drops of benzene forming in 0-52 N propionic 
acid). 
Taking ¢, (ie. & 1-67) for benzoic acid 
system, 
1 
a 


2.80 » < 36 
2-69 300) 0.854 (16) 


10° & 3600 0-492 (17) 


7 











Drop Volume (m. 
Os 


a? 











, 


r 


J 10 50 


me of drop formation and withdrawal ec 


Fig. 5. Variation of K,, with time of drop formation (for drops of benzene forming in 0-O1LL N benzoic acid). 
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whence 
aS 0-068 
Kun 0354 ~ 0-492 
giving A, = 0-34 1b mol hr ft® (Ib mol ft*) 
and 
O39 2-8 


0-37 Ib mol hr ft? (Ib mol ft*) 


Ap p 5DS 





giving A, 


p 


where A,,, is overall coefficient based on water phase 
concentration with benzoic acid as-solute. 

The experimental values for drops formed and 
KL, = 0-4 


collapsed in 1-67 sees (4, — lsec) are A, 


and kK, 


, of 0-68, which are about twice those given 


hor 00d wf drows se 


Fig. 6. 


(for 


Variation of Ay with time of drop formation 
drops of benzene [0-07 ml] forming in (52 .N 


propionic ac idl), 


by the Hicsie theory. The index for the experimental 
07 

t, 

difficult to give a satisfactory explanation for these 


work is whereas the theory gives ¢,°°. It is 
differences but the actual process of drop formation 
may be more complex than appears and eddies may 
give rise to more rapid transfer. 

2. Effect of drop size—The percentage equilibrium 
reached by the drop increases as the drop volume 
decreases. This is believed to occur because of the 
increased area of interface per unit volume with the 
smaller drops. For the range in drop volume ex- 
amined (0-04-0-114 ml), the overall transfer coefficient 
was found to be constant, but the range is too small 


for any definite conclusions to be reached. 


3. Effect of change of solute—The behaviour of both 


systems is basically similar, but with propionic acid 


as the solute the drops attain more closely to equi- 
librium with the water phase. The values of A for 
both systems are also similar. Since the resistance 
to transfer is in both phases with propionic acid 


as solute and mainly in the water film with benzoic 
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acid as solute, the latter system might be expected 
to give the higher value of A. On the other hand 
the diffusivity of propionic acid in benzene is larger 
than that of benzoic acid in water, and it is believed 
that it is this difference in diffusivities that causes 
the differences in the film transfer coefficient for the 
two systems. 


, 


Transfer to rising drops. Series 2 and 3 


Two series of experiments have been made in which 


the transfer to rising drops was measured. In series 2, 


& 
Verd/e 


vowe 


Fig. 7. Apparatus for measurement of extraction into 


rising drops, 


the drops were formed from a glass jet and rose 
through a column containing the continuous phase 
at rest. In series 3, arrangements were made for the 
continuous phase to flow down the column counter- 


current to the rising drops. 


Series 2. 

Transfer to single drops rising through still liquid 
The glass apparatus used is illustrated in Fig. 7. 
The column was | in. diameter and was closed at 
the bottom by a rubber bung through which passed 
the glass jet. The phase to be dispersed was contained 
in a large funnel and passed via a stainless steel 
needle valve through the jet into the column. A 
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range of tubes of various heights were available and 
several jets were used to give drops of different 
size. The dispersed drops on reaching the top of the 
column were directed through the inverted funnel out 
of the column and into a burette. In operation the 
column was filled with the aqueous phase to the neck 
of the top funnel and the needle valve was adjusted 
to give the desired rate of drop formation. It was 
important to maintain the top interface at a fixed 
position inside the funnel to give reproducible end 
conditions. 

The 


influence on the approach to equilibrium and on K 


following factors were varied and _ their 


were determined. 
1. Time of rise of the drop. 
Size of drop. 
3. Rate of drop formation. 
4. Physical properties and distribution relation 
the system. 
5. End effects. 


To obtain the value of A for the period of rise 
of the drops, experiments were made with two 
columns of 33 and 125 em high; by subtracting the 
results in the short tube from those in the long one, 
the end effects are largely eliminated, and A for the 
period of rise can be found. These experiments are 
conveniently arranged as: 

Series 2 A-E Two heights of column with five jet 
sizes, using Benzene-benzoic acid-water. 

F Three columns of differing diameter 

with the same jet size. 
Two heights of column with one jet 
using the system Benzene-propionic 
acid-water. 
As G but using Methylisobutylketone- 
benzoic acid-water. 


Calculation of transfer coefficient K for rising drops 
Using Benzoic acid as the solute the percentage 
equilibrium attained by the drop with the aqueous 
phase is very small, and taking the distribution 
coefficient as constant for any one test, A is given 
by, 
N 
S.ACm 


»» 


where N is the mols transferred in time ¢. It is, 
however, not so much the height of rise that is 
important but the actual time of rise. The method 


used can be set out as below for tests )) where the 
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time of drop formation was l sec and the drop 


volume 0-097 ml. (Interface = 1-015 sq. em.) 


Water phase conc. = 0-O11 N 





Time of rise Height of rise Benzene Equil. water 


normality normality 


aca cme 


12-05 125 


32 33 


O-O072 
0-0031 


O-0018 
O-OOTI 





O-OOTS — O-OOLT 


Hence 10, = O-O11 OMIIDS 


S600 
SSD 


2-5 « O- 162 » 097 
and = = 
LOLS © O-00955 


16-3 gm mols hr cm? (gm mols cc) 


0-535 Ib mol hr ft? (Ib mol ft*). 


The experimental results are presented in a series 


of tables and figures which are discussed below. 





| 





4 
| 
| 
| 

| 

| 


800 900 1000 





100 #200 300 YOO S00 600 700 
Reynolds Number 


Fig. 8. Variation of A, with Reynoips number (for 


drops of benzene rising through 0-011 N’ benzoic acid; 
height of rise 33-125 em). 


The influence of the time of rise of the drops is 
brought out together with the size of the drops 
since this latter factor alters the velocity of rise 
and K,, is determined by the difference in extraction 
over the two heights of 33 and 125 cms. 


Effect of drop size on K for rising drops. Series 2 A-E 
Fig. 8 shows the values of A,, plotted against Re 
based on the physical properties of the continuous 
phase and the drop diameter, and Fig. 9 values of 
K,, and the drop diameter. The averaged values for 
each series are given in Table 2 which shows the drop 
diameter and the velocity of rise. The values of K,, 
have been determined as explained and only refer 
to the period of rise. 
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Table 2. Experiments: Series 2 A-E 
Variation of A,, with drop size for rising drops. 
System Benzene-benzoic acid-water with O-O11 N  Benzoic 


acid as water phase. 





K w 
lb mols 
A r fe 


(lb mols/ft*) 


Drop V eloc wy 


“ luce 


Drop 


diameter ol res 


ml mm cm #ec 


O34 
(42 
O475 
535 


O45 


(horSs6 
O-Or6 
O4MS 
O97 
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Since an increase in the drop size does not of 


self alter any of the factors on which A’ is based, 





— 





i 
2 
Drop aameter 





Fig. %. Variation of A, with drop diameter (for drops 


of benzene rising through O-OLL N benzoic acid: height 


of rise 33-125 em). 


with Re be attributed to 


alteration of the flow conditions around the drop. 


the increase in K must 
The points on Fig. 8 lie on a straight line, whilst 
those of Fig. 9 lie on a curve. This arises from the 
fact that with very small drops the velocity of rise 
is less than with the larger drops and the Re group 
allows for this factor. For very low values of Re 
the drop should behave as a solid sphere and the 
On this basis 


flow around it should be streamline. 


a rapid falling off of the curve in Fig. 9 is expected 
values of Re. 


fi vr ke ~~ 


Ejfect of rate of drop formation on transfer 

The results of series 2 D are given in Table 3, from 
which it is seen that there is no significant effect 
of the rate of formation of the drop on the transfer 
occuring during rise. The total transfer will be affected 
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Table 3. 


Effect of time of drop formation on mass transfer to 


Experiments; Series 2D 


rising drops. 
System Benzene-benzoic acid-water, 


yet diameter 1-24 mm. 





Time of PF j e 
ned Wate r © 


normality 


Time of Drop 
drop 5 , Pe 
. ° rise volume he niene Equi - 
formation ad 

normality librium 


NOON SECS ml 

Column height 33 cm 
0-002 76 OO 
OOoll 
moll 
OLLI 
Oolll 
Olli 
OLLI 


3-2 O85 
3-2 O-OS | 00278 
3-2 0-092 


3-2 0-083 


0-00276 
0-024 

32 0-086 O-O0295 
3-2 0-098 0-00286 


3.9 O4V97 O-00306 


Column height 125 em 
Od] 


0-086 


O-OLIS 
OOLLS 
O-OLIS 
OOLLS 
O-OLLS 
OOLLS 
OOLLS 


0-062 
O-O0T6 
0-0069 
O-0073 


O-OS75 
0-004 

0-098 0-070 
O-O074 


O-097 
0-099 





since the transfer during drop formation is influenced 
by the time of drop formation. 


Effect of distribution coefficient and solute diffusiv- 
ity. 2G 

The the 
coefficients for the system benzene-propionic acid- 


method used for calculating extraction 
water was similar to that described above excepting 
that the drops approached more closely to equilibrium, 
the the 


different in the 35 cm column from that in the 125 em 


and extraction across top interface was 
column. The transfer was assumed proportional to 
the concentration at the top and the results were 
obtained by successive approximations. 

Results for the tests with 33 and 125 em columns 
are given in Table 4. For a drop volume of 0-072 ml 
with 


K,, — 0-5 for benzoic acid and Ay, = 0-525 


propionic acid. These values for rising drops are 
related to each other much in the same way as those 


given earlier for forming drops. 


Ejfect of interfacial tension 


To investigate the effect of a marked change in the 
interfacial tension the system methyl isobutyl-ketone 
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Table 4. Experiments: Series 2 Gand H 


Values of K,, for rising drops. 
Tests G. System Benzene-propionic acid-water. 


Tests H. System Methyl isobutyl ketone benzoic acid-water. 





K 


qm mols 


w 


Drop 


Tests volume diameter 


Drop Velocity 


lb molhr iC 
lb mol /j® 


of rise 
hr cm* 


ml mm cm sec qm mol ce 


G oo72 5-16 
H 0046 4-4 


10-4 560 16-1 
95 440) 49-4 


0-53 
1-62 





In this case EZ is 
very small and A,, is calculated as on p. The ex- 


With both 


the systems methyl isobutyl ketone benzoic acid. 


benzoic acid-water was studied. 
perimental results are given on Table 4. 


water and benzene-benzoic acid-water the distribution 
coefficient is very low and most of the resistance 
will lie in the water phase. One would 
K,. The 


values for A, with methyl ‘sobutyl-ketone as the 


therefore 
expect similar values for experimental 
solvent are, however, three times those with benzene. 
These differences cannot be explained on the existing 
simple ideas of diffusion and other factors must be 
involved. The interfacial tension between M.1I.K. and 
water is only &Sdynes em while for benzene-water 
it is 35 dynes cm. This is a large difference and may 
be the cause of the big difference in A,, for the two 
systems. Further evidence on the influence of inter- 


facial tension will be given in Part 2. 


Studies of the wall effect. Series 2F 

In the experiments the columns used were of suffi- 
ciently small diameter relative to the drop to influence 
to some extent the velocity of rise of the drop. It 
the effect of 
tube size in relation to the drop diameter. The results 


was therefore decided to investigate 


as given in Table 5 show clearly that A is  inde- 


Table 5. 
Studies of the wall effect. 


Experiments; Series 2F 


Svstem Benzene-benzoic acid-water. 


Drop diameter 5-5 mm. 





Velocity of rise 7 
lb mol/hr fC 
lbh mol /t® 


Column Diameter 


cm in. 
cm sec 


ft hr 
O39 
O50 
1-18 
1-18 


6-74 805 
8-25 O75 
10-6 1250 
10-4 1230 


O-505 
0-49 
47 


0-52 





K,, calculated using data for end effects given in Table 6. 
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pendent of the tube size for drops of the same diameter. 


The drop as it rises through its own height displaces 
an equal volume of the continuous phase between 
itself and the tube wall and the relative velocity 
between the drop and the continuous phase is un- 


altered, although its absolute velocity is reduced. 


Calculation of end effects 
Some indication of the end effects can be obtained 
by calculating the anticipated extraction for a rise 
of 0-33 cm in 3-2 secs, based on KA, obtained for 
the rising section only, and comparing this value 
with the experimental figure. 

For a period of rise of 8-85 sees (p. 205) 


E oe). 


Hence for a period of rise of 32 SsCCS, the expected 


value of E is D5 R9 


ODL, 

The experimental figure for this transfer is 2 
1-725, so that the end effects give rise to a value of 
O-S1E£, or nearly 47% of the total transfer occuring 
in the column of 33cm. If corrections are made for 
the small change in driving force between the section 
still be shown that 

10% of the total 
transfer for the 33cm column. A summary of these 


0-33 em and 33-125 em, it can 


the end effects amount to about 


results is given in Table 6. 


Table 6. Experiments: Series 2C EK 


Summary of end effects for mass transfer to rising drops. 


System Benzene-benzoic acid-water. 





Time of drop formation 1 sex 
Series 2C: Jet 0-78 mm 
0-049 


0-405 D2 


Drop volume mi. O-0495 


“ Equilibrium for formation 
% Equilibrium for both ends. 1-0 
DU’ »\ ee: as a ss os Dee 
1-24 mm 


Drop volume . 


‘eries 2D: Jet 
O55  O-OS 


*» Equilibrium for formation O33 0-38 


% Equilibrium for both ends. 0-76 Oso 
Factor R 
‘cries 2 E: Jet 3-65 mm 
Drop volume ml. 
% Equilibrium for formation 


‘» Equilibrium for both ends . 
Factor R 





R Combined end effects 
Transfer during drop formation — 
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Nerves 3. 
Extraction into drops rising through a moving liquor 


In a number of practical examples of extraction the 


drops rise against a falling stream of the continuous 


phase. These conditions have been produced using 


© 
© Aipoedie > 
Y 


7ve 


Fig. 10. mass transfer into 


drops rising through 


Apparatus for measuring 


a moving contimuous phase. 


The 


continuous phase was controlled by 


flow of the 


the 


the apparatus shown in Fig. 10. 
stainless 
steel valve and the interface maintained in the neck 


of the funnel at the top. The drops were collected 
Table 7. 


Exype riments: Neriea 3 


Effect of velocity of continuous phase on transfer 


System Benzene-benzoic acid-water. 


Jet diameter 0-114 mm. 


Height of moving column SO cms. 
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in a small reservoir above the column. In operation 
the continuous phase was first started, and when it 
was flowing at the desired rate the dispersed phase 
was introduced. The experiments were made with 
the system benzene-benzoic acid-water, and the time 
of drop formation kept constant at 2 secs. 

The results are presented in Table 7, which gives 


the mean downward velocity of the continuous 
phase, the absolute velocity of rise of the drops and 
the total time of rise. Under these conditions, Re 
based on the relative velocity of the two phases and 


the drop diameter does not alter, any increase in 


e of e of & 


Fig. 11. 


(for drops of benzene rising through moving liqnaid). 


Graph of “o equilibrium vs. time of contact 


the downward velocity of the continuous phases 
being equivalent to a decrease in the upward velocity 
of the drop. The important effect of an increase in 
the downward velocity of the continuous phase is 
that it 
of contact of 


increases the time 
the drop 
to rising drops. with the liquid for a given 
Thus the 


results are conveniently ex- 


height of rise. 


pressed for this series by 





Time of 
Water 


” loe ity 


Time 


of rise 


Drop 


volume 


ly op 
drop 
1% lo« ily 
formation 


arcs cm acc cm ace 


th (-Ors4 
moo 
(F085 
101 
mal 
104 
LIS 
116 


O-OLLG 
O-OLL7 
O-O109 
O-O108 
O-OL09 
OoLl2 
O-OLLS 
0-0106 


plotting the approach to 
the 


Final 


he niene 


Water 
nor mality 


equilibrium against 


normality time of contact, as shown 
in Fig. 11. This graph in- 
0-0064 dicates that FE varies with 
0-007 
0-074 
0-008 1 
0-0077 
0-0092 
oll 
0-0086 


the time of contact of the 
two phases. 

Fig. 11 also brings out 
the fact that 
transfer is increased equally 


the actual 


by increasing the height of 
the column with zero flow 
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of the continuous phase or by increasing the time of 
contact to the same extent by increasing the down- 
ward velocity of the continuous phase. Since AC), 
and the drop volume are approximately constant the 
linear relationship between EF and the time of contact 
that the transfer coefficient is in- 


is an indication 


dependent of the latter. 


Mechanism of mass transfer to rising drops 


The results of the experimental work show that: 
(1) A is independent of the time of contact. 
(2) A increase with increasing drop size. 

(3) A increases with decreasing interfacial tension. 

(4) For the systems tested, A for a rising drop 
corresponds with A for a drop forming in a given 
time, irrespective of solute. 

Consider the interface of the rising drop. At any 


time, ¢, a point, A, on the interface of the drop, a 


point, B, in the surrounding medium, and (, the 
centre of the drop are in the same horizontal plane. 
At a time, f 
for point 


At, there are three possible positions 

A in relation to B and ©. 

(1) A and (' remain in the same position relative to 
each other, the drop behaving as a solid sphere. 

(2) A remains in the same position relative to B (Ag) 

For this to 


move downwards relative to the drop centre 


in the diagram. occur, A must 
at a velocity equal to the velocity of rise of 
the drop. The drop interface is therefore being 
continually formed at the top of the drop and 
destroyed at the bottom. Thus, the interface 
will be replaced every time the drop rises 
through its own height. 
moves to some point, 4g, between A, and Ag. 
The drop is continually renewing its interface, 
but the tangential velocity at the equator of 
the drop is less than the velocity of rise. 
Case (1) was considered by SHERwoop [7], who 
showed that the predicted extraction was much 
If we take the 


case of propionic acid, where the water phase resist- 


smaller than the practical results. 


ance may be neglected, then the anticipated extrac- 
tion may be calculated by NewMan’s method [12]. 
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Thus, 


time of formation and rise (sees): 
radius of drop: 
mean concentration within drop: 


an integer number. 


For a time of drop formation and rise of 14 sees 


(Tests @) Dy 1-72 cm? sec. 


DO 1-72 4™« 10° 


—— 0-356 . 
as 2 


Drop diameter == 0-52. em (Table 4). 


From NEWMAN'S paper, the corresponding value of 
Ww 
n 


‘ 


is OSI. This gives the predicted equilibrium 
as 19%, whereas the experimental figure is 55%, or 
nearly three times as much. 
the 


If the water phase 


resistance were included difference would be 
vreater. 
In examining case (2), we may apply Hicsie's 


equat ion 


for a velocity of rise of 10-4 cm sec 


Diameter of drop — 0-52 cm. 


ky, 
10°? gm mols hr em*® (gm mols cc) 
2.5 lb mols hr ft? (!b mols ft*) 


OSS & 1LO-% 
OOD 


k, 113 | 


1-5 


1-78 lb mols hr ft? (Ib mols ft*) 


10°? gm mols hr em? gm mols ce) 


whence 


Ap 2-76 . 1-78 


K ,, = 1-66 lb mols hr ft® (Ib mols ft*). 


This result is three times the experimental figure, 
which lies between case (1) and (2). Thus, the mechan- 
ism may be considered as represented by case (3), 
in which the drop renews its interface as it rises, 
but that the tangential velocity of the interface is 
less than the velocity of rise. The rate of interface 
formation per unit of interface would be expected to 
decrease with drop size and increase with decreasing 
interfacial tension. The general results are in agree- 
ment with these two features. 
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Application of data for single drops to spray and 
perforated plate towers 

The experiments so far described have provided data 
on the extraction occuring during formation, rise and 
coalescence of the drops as the top interface. These 
three stages form the basis of the operation of spray 
and perforated plate towers, and it is therefore of 
interest to examine the performance of such equip- 
ment in the light of the data obtained. 

ALLerTON, Strom and Treysat [1], have pub- 
lished their results on the system toluene-benzoic acid 
water using a perforated plate column of the foll- 
owing construction. 

Overall tower height . . 5 ft 3) in. 
Internal diameter 3°63 in. 
Cross sectional area . 0719 ft? or 66-6 em*, 
Plates 11. Thickness 1 16 in. Spacing 4-75 in. Holes 3:16 in, 
diameter. No. per plate 51. Disengaging space at top 7-625 in. 

In operation the extraction was from toluene to 
the water, the toluene entering at the bottom. Since 
the distribution coefficient for benzoic acid in benzene 
water is 0-068 and for toluene water is 0-070 and the 
physical properties of the two systems are very 
similar the present data for benzene drops will be 
applied to examine the toluene system in the plate 
column described above. 

Total 


SOSTS in 


height of column 75.11 


Effective volume of column MOT1o 
O-358 ft*, 
Flow 


32 cc sec, 


rates: Toluene 57 ft® ft® br 4-1 ft® he 


Ti . St. Deop volun 

Time of drop formation 4 

SL OST 
32 


Water 75-6 ft® ft® hr 5-44 ft® br 


The transfer coefficient for the three stages of forma- 


0-50 secs *. 


42-6) ce sec. 


tion rise, and coalescence can now be calculated and 
hence the overall transfer coefficient deduced and com- 
pared with the experimental value given by ALLERTON. 
Forming drops—For a drop of 0-37 ml forming in 0-59 secs K,, 
can be obtained from Fig. 5 as 1-3 1b mola hr ft*. 

1-49 om? 

149. 51 

3.5? 

O-Om2 ft? 


Average interfacial area of 1 drop 


Total area of forming drops on one plate 


Total area of forming drops in columns — 082 » 11 

(re tte 

(K . 8), forming drops 3) O-oR 
1-17. 

* Obtained from data for drop volume vs.time for forma- 


tion for ed holes. 
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Rising drops—Allow as top interface at each plate a layer 
of jin. as given by authors. 

13 cm see 

0-64 cm sec 
12-4 cm ‘sec 
2-46 secs. 


+ 7-625 = 47-62 in. 


Velocity of drop in still water . 
Velocity of water phase 42-6 66-6 
Actual velocity of rise of drop 
Time of rise for 1 ft ° 
But height of column of rising drops 10 » 
47-62 « 2-46 
12 


845. 


Total time of rise 


No. of drops oS x 5l 


: 
5a 


S45 2-40 
30-52 


for rising drops of 0-37 ml at 12-4 em/see can be 


Hence total area of rising drops S ;, 


But K,, 
read from Fig. 9 as 0-69 lb mols hr ft®. 


Hence (AK. 8), rising drops — 2:26 « 0-69 — 156, 


Top interjace—For drops of this size the tranfer occuring 
at the top of each plate may be taken as equal to that in 
formation (see Table 6). 

The 


whole tower is then given by: 


overall transfer coefficient as measured for the 


K..a.V 
_# .a ~ 0358 


K,, a 


1-17 1-56 1-17 30 
10-0 Ib mol hr ft® (Ib mol ft*). 


Similar calculations have been made for two 


other flow rates and the results together with the 
experimental value of TREYBALL are given in the 
table, 





Water tlow K 
eiehr 


Toluene flou 


ee hr 


ae@ Ky, .@ 


calculated experimental 





Similar calculations have been made on the data 
of Treysat and Dumovu iy [9], using a perforated 
plate column very similar to that described above 
but in which the plate spacing could be varied. For 
three different flow rates, the values of A,a were 
12-6, 7-4 and 7-4, whilst the calculated value based 
on this research was 12-7, 8-7, and 82. These cal- 
culated values are in very good agreement with the 
experimental values and it seems reasonable to 
conclude that the columns do in fact behave much 


as indicated by this work on individual drops. 


SUMMARY 
The extraction occuring from a continuous phase to 
dispersed drops has been measured both for condi- 
tions where the continuous phase is stationary and 
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when it is flowing countercurrently to the drops. 
The three stages of transfer have been measured 
directly, that during drop formation, that during 
the rise of the drop, and finally that occurring at the 
top interface of a column where the drops coalesce. 

Mass transfer to a drop during formation has 
been found to be almost independent of the time of 
formation of the drop for a range of from 4 to 1 sec 
for formation. The overall transfer coefficient K 
based on the average area exposed during formation 
of the drop decreases with increasing time of forma- 
tion, but is practically independent of the drop size. 
Smaller drops approach more closely to equilibrium 
because of the increased area of interface per unit 
volume. 

K for transfer during the rise of the drops has 
been found to increase with drop size, decrease with 
increasing interfacial tension, and is directly in- 
fluenced by the relative velocity of the drops and the 
continuous phase. It is the time of contact rather 
than the height of rise that is important. Comparison 
of the experimental work with the anticipated results 
assuming the drops to rise as solid spheres or alter- 
natively that they renew their surface each time 
they rise through one diameter suggests that the 
It is, 
however, felt that the reformation of the interface 


real conditions lie between these two cases. 


together with any recirculation of the liquor inside 
or around the drop that is the determining feature 
in the process. For a further understanding of the 
mechanism of the process measurements on the 
velocity of the interface are required. This is very 
difficult with drops but in Part 2 of this work the 
effect of the velocity of the interface in a horizontal 
extractor with the two streams flowing side by side 
in a tube will be discussed. 
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NoTaTION 
a = Area of interface per unit volume of apparatus (eq. 5) 
C = Concentration of diffusing solute (mols/ft*) 
C, ~ Concentration of diffusing component in phase A 
(mols/ft*) 
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C,, ~ Concentration of diffusing component in phase B 
(mols/ft*) 

Concentration of diffusing component at the inter- 
face in phases A and B 


Cir Ovi 
Concentration of diffusing component in phase 
in equilibrium with concentration C, in phase 
Concentration of diffusing component in phase 
in equilibrium with concentration C, in phase 
Concentration of absorbed gas in liquid phase under 
equilibrium conditions 
Initial concentration of gas in liquid phase 
Logarithmic driving force in concentration units 
Diffusion coefficient (normally in ft*/hr) 
Percentage approach to equilibrium 
Ratio of concentration of solute in equilibrium in 
both phases 
Film transfer coefficients. mols/hr. unit area. (unit 
concentration difference) 

Film transfer coefficient based on concentration in 

the aqueous phase 

Overall transfer coefficient |b mols hr ft? (Ib mols/ft* 

Overall transfer coefficient based on aqueous phase 

concentrations Ib mols/hr ft? (lb mols per ft*) 

Overall transfer coefficient based on concentrations 
Subscripts 
w and B denote water and benzene respectively. 


in phase A and phase B respectively. 


Mols of diffusing component transferred 

Time (normally in hrs) 

Time of exposure of interface 

Time of formation and collapse of drop 

Time of one cycle for liquid to enter and leave 
drop apparatus 

Radius of drop 

Area of interface 

Volume of liquor fed to jet that does not enter 
the drop 


Volume of drop 
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Summary 


The conditions of heat transfer to a fusible solid in contact with a heated surface are examined 


A theoretical equation for the heat transfer coefficient is deduced, and compared with practical results. 


A shorter, empirical formula is suggested for design purposes. 


Résumé 


Etude de la transmission de la chaleur entre une surface chaude et un solide fusible 


Etablissement théorique d'une formule donnant le coefficient de transmission et comparaison aver les 


resultats experimentaux 


L'auteur recommande Pemploi d'une formule empirique simplifiée, suffisante pour la pratique indu 


strielle. 


Although the problem of heat transfer to fusible 
solids arises in process design, the literature contains 
few references to the subject. Mathematical solutions of 
special cases have been provided by several authors [1], 
but practical investigations have apparently been lim 
ited to instances of ice formation [2], or of solidification 
and fusion on immersed surfaces [3]. 

A fusible solid may be melted by one of the follow 
ing methods: 

(a) The solid may be added to a mass of already 
molten material in a heated, agitated vessel 

(b) The solid may be arranged so that hot gases 
(or its own superheated vapour), may be circulated 
over it, the molten portion draining away. 

(c) The solid may be placed on a heated surface, 
and the molten material drained away. 

In case (a) once enough liquid has been accumulated 
in the vessel, the heat transfer is via this liquid 
Hixson and Baum [4] have shown in their studies of 
mixing how this problem may be solved by correlation 
of the agitation, and physical properties of the system 

In (b) the heat necessary to raise the solid to melt 
ing point, and to bring about the fusion must be trans 
ferred through the gas film and liquid film successively 
Work on this mechanism is in progress. 

In example (c), the sensible and latent heats must 
be transferred from the hot surface to the melting 
solid through the film of molten matter. It is this 
problem which is examined here. 


The liquid film 

The condition of a film of liquid between two plane 
surfaces has been investigated, for example by Mrr- 
CHELL [5]. The present case is somewhat complicated 
by the fact that the fluid is not introduced along any 
preferred direction, but is generated at the solid surface. 
In general, however, the conditions of flow in such 
a film will be controlled by the viscosity and density 


of the liquid, the inclination of the film to the vertical, 
and the pressure exerted on the film by the weight of 
melting solid. 

Then, the difference between the frictional forces 
at the boundaries of a small layer of the film (of thick 
ness dx, and distant + from one face) must balance the 
accelerative forces when the flow rate is constant 
Thus 

dV div +dV) 


uM. “”. q\a 


j F 0 cos @) dx 
az . ar > 


for unit breadth and so, 
g(a 0 cos O) dx. 


If there is no slip at either the heated surface or that 
of the melting solid, and if the latter is also stationary, 


integration gives 


dx 
0 cos O)\m . ~ 


dV =—qle x.dx) 


where ““m” is the total film thickness. 

Substituting the mass flow rate at any point, 
(dW dt), for the velocity term, and integrating over 
the film results in the expression 

l2y dW 
go(a+ecos)*\ dt | 


(1) 


Heat transfer through the film 


The flow of heat through the film covering a small 
portion of the heating surface, dA, may be written as 


dQ : 
tt aA, T,,) (2) 


so that 
d Y dQ 


h dt.dA.(T; — Tm) dt. B.dL.(T, — Ty)" 


when the surface has constant width B. If the liquid 


flow rate at this point is dW dt also 


dQ aw 


Y= A+ 0,(T. — 7] (3) 
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Combining (2) and (3) 
dW[A+Cp(Tm—T)] 
dt. B.dL (Ts — Ty) 
k 
m 

if the flow is viscous. 

This expression is comparable to that developed by 
Nussevt [6] for condensate films, and a similar process 
of integration over the entire surface results in the 
following expression for the mean film coefficient, h,, . 


kM gola- 0 cos 0) | 7 . Ce (F 
h 0-666 = : i 
oa a. LT, T..) 


or, in dimensionless form, 


0 00s 9) [A+ Cp (Tm 
kp. (TT, — Tn) 


Lgqola , T 025 


he. L 
h 


0-666 . (6) 


If the heating surface is provided as a horizontal tube, 


of outer diameter D,, L becomes and @ has 


a minimum value of 0°, so the alternative formula is 


then 


= LE gola- 2 OLA+Cp(T, —T))}}025 2 
0-710| 2o9el + geos OA + Cy P97) 


hm. Do 
keg. (7, — Tm) 


k 
Now eqs. (5), (6) and (7) apply only when the flow in 
the film is viscous, and the conditions for turbulence 
will now be examined. For turbulence to develop, on 
a plane surface, 

,.dW 2 F 
Re ~ 2100 - 4mB.dW (<= ). 
2(m+ B).dt.mBu Bu \ dt 
when m is small, 


dw 


must exceed 
dt 


or 1050 Bu. (8) 


The flow rate at the end of a length of heating path L. 
may be written as (dW/dt), where, 


dw hn (Ts Tm) BL 


9 
| dt L [A T Cp (Tm T)| ; \ 


Thus, for turbulence to occur, at this point, by combin- 
ing (8) and (9) 


hm (Ts — m) BL 


: O50 
[A+ Cp (Tm —T)) >I Bu 


. 1050 il [A » Cp Tw» - T) 
L must exceed hn A(T, — Tn) ; 
Substituting for h,, from (5), 


L must exceed 
nla Op (Tm — F))})| w 
k(T; — Tm) | 


(*Substitute 23360 when D, replaces L). 


0-33 (11) 
18,350* 


ge (a + gcosé | 
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In all ordinary cases, therefore, the length Z is so small 
comparatively, as to preclude turbulence. 


Apparatus 

The apparatus consisted of a heating surface—either 
in the form of a steel platen, or copper tubes—arranged 
ina sump. The dimensions of the plane surface were 
varied between 10 em and 50cm in both length and 
breadth. Tubular surfaces were provided as two par- 


allel rows of tubes, arranged on a 60 degree pitch of 





ef 














Ow 08 } 16 20 
Log L 
Effect of length of heating surface upon heat 
transfer. 


Fig. 1. 
one outer tube diameter. Tubes with outer diameter 
of 0-635 em, 1-270 em, 2:154.em, and 5-080 cem were 
used. The heating surface could be fixed at any angle 
to the vertical by means of adjustable legs. Insulated 
side plates could be adjusted to enclose any desired 
amount of heating surface. The heating surface was 
fitted with thermocouple junctions. 


Procedure 
In an experiment, the heating surface was adjusted to 
the required angle, with the front edge horizontal. 
A charge of material, prepared as a solid block or as 
sized pellets, was added to the apparatus, thermo- 
couples inserted, and weights placed upon it, if it 
was desired to alter the loading. 

Steam was then admitted to the heating surface, 
and after an initial steadying period, the molten sub- 
stance was collected and weighed at timed intervals. 
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Fig. 2. Effect of loading pressure (¢) and inclination of 


-~. 


heating surface (@) upon heat transfer. 
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Fig 3. Effect of 7, — 7, upon heat transfer. 
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Fig. 4. Correlation of general results, by means of eq. (6) and (7). 


Results 
Mean heat transfer coefficients were calculated from 
the observed melting rate, the latent and sensible 
heats involved, and the temperature difference between 
the heating surface and the solid. 

The physical properties of the molten film were 
extracted from usual sources [7]. Much of this data 


9 


was not available at all temperatures, and extrapola- 
tion was used. 


Eqs. (6) and (7) suggest that the effect of the length 
of the melting path (LZ or D,), the bearing pressure 
and the draining effects (o + 0 cos 6), and the tem- 
perature difference (7; — 7),) should be investigated. 
Some results are presented in Figs. 1, 2 and 3, in which 
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the full lines represent the slope predicted by eqs. (6) 
and (7). 

These, and other results are plotted in Fig. 4 to 
show their relationship to the theoretical eqs. (6) 
or (7),—represented by the broken line. 

Although the scatter of individual points is serious, 
it will be observed that for any single compound the 
divergence is less wide. It would appear that the 
unreliable physical data is responsible for this. It is 
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or, 
hin. Do o cos 6) [7 T 
ik 


kw Ts 


4g ole T Cp (Tm " T)| ° ” 
| 


Tm) 


0-710 | t 


for tubes. 

The following shorter form of this equation was 
found to be satisfactory for use with tubes, and gran- 
ular organic materials. 


300 io aa 
h (7, — T,)°*8. De J.sec™.m=* . (°C). 














NOTATION 


heating area m* 

width of heating surface m 
Specific heat of fusible solid 
J kg. (°C) 


outer diameter of tube m 





acceleration due to gravity m.sec™? 
local heat transfer coefficient in molten 
film J.sec™!. m-?.(°C)" 

mean heat transfer coefficient in molten 


9 


film J.sec™!.m-?.(°C)"! 





Fig. 5. Comparison of some practical results with the empirical equation. 


"= (0, Ta D5 oe ere 
also of interest to note the degree of agreement be- 
tween results obtained on tubes, and plane surfaces 
It should also be recorded that when the solid material 
was prepared as pellets, the loading effect was appa- 
rently limited to that due to a height of charge of 
about three times the container side. 


The results for granular organic compounds, and 
tubular surfaces may also be interpreted by a shorter, 
empirical equation, in which those physical properties 
which are reasonably similar for different materials 
are neglected. This formula is 

300) 


= A 2 y 1 
p?™ (7, rs Tn)*® pers J . BEC ~ mn ~ ( ) ° 


(12) 
The correlation given by this equation may be judged 
from Fig. 5. 


Conclusions 


The theoretical equation for heat transfer to a fusible 
solid may be evaluated from the formulae for plane 
surfaces, 


L wo | LAgola + ocos6)|A+ Cpi Tm — T)| 0% 
= 0-666 = ; 
A k. (Ts — Tm) 


thermal conductivity of molten film 
J sec. m-!. (°C) 

length of plane heating surface m 
thickness of molten film m 

heat J 

REYNOLDs number 

time sec 

initial temperature of solid material “C 

melting point of solid material °C 

Temperature of heating surface °C 

velocity in film m.sec™ 

mass flow in film kg 

dimension in molten film m 

angle of heating surface to vertical 

latent heat of fusion J.kg™ 

viscosity of molten film kg.m-'. sec! 

density of molten film kg .m~* 

pressure gradient in film due to weight of solid 
kg .m-* 
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Résumé 


L’auteur donne un apergu des problémes technologiques que pose la mise au point des procédés 


catalytiques, suivant les divers modes operatoires choisis: travail sur lit fixe, travail par percolation, 


travail sur bouillie, travail sur lit mobile ou sur lit fluidisé. 


Les résultats des études fondamentales montrent que celles qui sont relatives au lit fixe sont les plus 


avancées, mais qu'il nous manque encore beaucoup de données pour les autres modes opératoires. 


Summary 


A survey is given of some chemical engineering problems occurring in the development and 


design of catalytic processes. They refer to the various methods for carrying out catalytic reactions, viz 


fixed bed, trickle, slurry, moving bed and fluidized bed operation. 


The results of fundamental studies show, that so far most data are available for the fixed bed processes 


but that many more are required particularly when referring to the other methods of operation. 


1. INrropvcTion 

Le développement des procédés catalytiques dans 
Vindustrie du pétrole et des carburants synthétiques 
est trés remarquable. C'est aux Etats Unis surtout 
qu'on s'est rendu compte de l'importance de la cata- 
lvse hétérogéne et particuliérement de la catalyse de 
contact. Ce développement est d’autant plus frappant 
si nous considérons l’accroissement des procédés pour 
la fabrication d’essence d’aviation. Par exemple, pour 
le cracking catalytique, la capacité est passée de 10000 
barils par jour en 1936 a 1700000 barils par jour en 
1950, 

La solution des problémes qui se sont présentés et 
qui se présentent toujours dans la réalisation d'un 
procédé catalytique est seulement possible s'il y a une 
collaboration incessante entre le chimiste, faisant les 
recherches fondamentales concernant le catalyseur, et 
Vingénieur du génie chimique. 

En ce qui concerne les qualités exigées d'un cata- 
lvseur industriel, celles-ci ont été décrites largement 
par Hooe [8] et peuvent étre résumées comme suit: 

1) Activité: les conditions de réaction doivent 
étre telles que léquilibre soit le plus favorable. 

2) Sélectivité: la composition du catalyseur doit 
étre telle que la réaction qui nous donne le produit 
désiré soit favorisée. 

3) Durée de vie du catalyseur: ceci est de 
premiére importance pour l'économie du_ procédé. 

4) Résistance mécanique: ce facteur a une im- 
portance spéciale pour les convertisseurs A lit fixe. 

Dans cette communication je me propose de vous 
soumettre un apercu d'un certain nombre de problémes 


de technologie dans la mise au point d’un procédé 
catalytique, aussi bien dans lusine pilote qu’a l’échelle 


industrielle. 


La NATURE DES PROBLEMES ET LES DIFFERENTS 
Mopes OPERATOIRES 


La nature des problémes qui peuvent se présenter 
a conduit A différents modes operatoires pour les 
réactions catalytiques. La conception des convert isseurs 
demande en général une connaissance des facteurs 
suivants, qu'il ne faut surtout pas sous-estimer puis- 
quils déterminent les possi bilités du procédé au point 


de vue économique, (Ce sont 


1) les lois hydrodynamiques qui portent sur 
écoulement des fluides et du catalyseur; 

2) les lois de transfert de matiére entre les 
fluides et les particules de catalyseur 

3) les lois de transfert de chaleur du catalyseur 


aux parois du réacteur et vice versa. 


A cété de ceux-ci on peut avoir le probléme de 
la régénération du catalyseur, par exemple dans le 
cas du dépét sur le catalyseur d’un produit dimi- 
nuant son activité. Cet effet est trés marqué dans 
le cracking catalytique, ob le dépot de coke peut étre 
de Vordre de 7 gl. La régénération se fait par brailage 
du coke avec de lair éventuellement dilué avec de 
lazote ou de la vapeur d'eau. 

Les différentes méthodes pour réaliser des réactions 
catalytiques a échelle industrielle sont résumées dans 
le tableau suivant et illustrées par un exemple in- 
dustriel. 
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Tableau 1. Modes opératoires des procédés catalytiques 





Travail sur lit fixe Travail sur touillic Travail sur lit mobile Travail sur lit fluidisé 


A. Réaction en phase B. Travail par percolation | Systéme: gaz-liquide Réaction en phase Réaction en phase 


yazeuse Systéme: gaz-liquide Catalyseur en suspension gazeuse gazeuse 


Liquide sous forme de dans le liquide 


Couche mince 


E xemples Exemples E xemples Exemples Exemples 


Hydrogénation des Hydrogénation des 


oléfines oléfines 


Hydrogénation du 


charbon 


Oxydation de Oxydation de 


l'o-xyléne l’o-xyléne 
Oxydation de 


paraffines 


Cracking du 
gasoil 
Synthése des 


hydrocarbures 


Synthése des 


hvdrocarbures 


Synthese des 


hydrocarbures 


Cracking du 
gasoil 
Synthése des 


hvdrocarbures 


Cracking du 
gasoil 
Synthése des 


hydrocarbures 





Nous ne voulons signaler ici que quelques procédés 
qui peuvent étre réalisés de facgons différentes. Les 
exemples comprennent des réactions d hydrogénation, 
d’oxydation, de cracking et de synthése des hydrocar- 
bures. 

Dans la section suivante on va examiner les pro- 
blémes qui se présentent pour les différents modes opé- 
ratoires. Quoique les développements récents nous 
aient donné les procédés sur lit fluidisé et sur lit 
mobile, c'est surtout le réacteur sur lit fixe qui est 
toujours le plus utilisé A cause de sa simplicité. 


3. Travam sur Lit Fixe 
A. Réactions en phase gazeuse 


La conception d’un réacteur sur lit fixe n’offre aucun 


probléme sérieux si la réaction a lieu avec un effet 


thermique peu important. Cependant pour la plupart 
des réactions catalytiques, effet thermique est consi- 
dérable et élaboration du réacteur le mieux adapté 
demande une étude importante. 

Nous n’avons qu’a songer au cracking Houdry, ot 
pendant une période de régénération de 10 minutes 
dans un réacteur ayant un volume net de catalyseur 
de 30m, la quantité de chaleur développée atteint 
1.3 10° keal. C’est surtout ce probléme de transfert 
de quantités de chaleur importantes qui a conduit au 
développement du catalyseur fluidisé. 


Nous pouvons maintenant distinguer les problémes 
qui ont trait a écoulement des fluides, et ceux con- 
cernant la transmission de chaleur et le transfert de 


matiére. 


a) La perte de charge dans un lit fixe—La perte 
de charge pour l’écoulement d’un fluide A travers un 
lit fixe de catalyseur peut étre caleulée au moyen de la 
relation existant entre la perte de charge et le facteur 
de frottement, celui-ci étant une fonction du nombre 
de ReyNnoups. Cette relation est analogue A celle pour 
écoulement d’un fluide dans les conduites. D’aprés 
BLAKE [3] et CaRMAN [5] la relation est la suivante: 
pet gp = 1(Re) 
ou 
lp — perte de charge 
o — densité du fluide 
v — vélocité linéaire 
1 = hauteur du lit 
S = surface extérieure du catalyseur par unité 
de volume du réacteur 


€ = porosité du lit. 
Un apergu des données pour des catalyseurs de 


forme différente a été donné par VERscHOOR [14]. 


b) Transfert de chaleur et de matiére—Ce transfert 
se fait par conduction thermique respectivement par 
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diffusion et par convection. Suivant les auteurs, 
c'est tantét la transmission de chaleur par con- 
duction, tantét le transfert par convection qui est 
le plus important dans la transmission de la chaleur 
aux parois du convertisseur A lit fixe (transport) 
normal au sens de l’écoulement du fluide). C'est grace 
au travail de Bakuurov et Boreskov [1] qu'on a 
appris que le transfert a lieu aussi bien par conduction 
que par convection, DAMKOHLER [7] a introduit pour 
le transfert de chaleur et de matiére dans un lit fixe 
les concepts de conductivité tourbillonnaire (eddy 
conductivity) et de diffusivité tourbillonnaire (eddy 
diffusivity), en relation directe avec la vitesse linéaire 
du fluide v ainsi qu’avec le diamétre du catalyseur 
d, suivant: 


sans dimension. FE a la 
(L2 T°). 


Pour le cas analogue de la conduction thermique, EF 


ou C est une constante 


dimension d’un coefficient de diffusion 
doit étre multiplié par c,o, oi ¢, est la chaleur 
spécifique et o la densité. DamKOnLeR donne la 
relation suivante: 


hotter til Ay Cp oO ( ; v d,, 


ou A, est la conduction thermique du lit sans écoule- 
ment de fluide. 

C'est A KLINKENBERG [10} que nous devons d’avoir 
montré, en dérivant la valeur de C des différentes 
recherches expérimentales, qu'il y a une bonne corréla- 
tion entre la théorie de DAMKOHLER et les données de 
Bakuvurov et Boreskov, de Bernarp et WiLneLo [2 
et de Coperty et Marsnauy [6): 





Vale ur 
de C 


Auteurs Transtert de 


Fluide 


Bakuvuroyv et Bo- 


RESKOV Chaleur et matiére Gaz 


Berxarp et Wu 
HELM Mati¢re Liquide 
COBERLY et 
MARSHALL 


Chaleur Gaz 





B. Travail par percolation 


Dans ce cas, un film d'un liquide traverse le lit de 
catalyseur, le gaz passant généralement en courant 
paralléle. Les problémes, par exemple la distribu- 
tion du liquide et la teneur en liquide (le « hold-up») 
de la colonne, sont surtout de nature hydrodynamique 
et en général analogues A ceux qu'on rencontre dans 
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les colonnes de distillation ou d'absorption avec an- 
neaux Rascuic. De cette fagon on a trouvé que la 
distribution uniforme du liquide dans le lit de cataly- 
seur s'effectue seulement si le rapport entre le dia- 
métre du réacteur et celui de la particule est superieur 
& une valeur minimum. 


4. Travait sur Boviiur 

Dans ce cas, le gaz traverse en forme de bulles 
une couche de liquide dans laquelle le catalyseur 
pulvérulent se trouve en suspension. L’application de 
cette opération pose surtout des problémes hydro- 
dynamiques comme la distribution des bulles de gaz 
dans la colonne de liquide et la formation des bulles 
ainsi que leur grandeur. Une étude préliminaire de 
Verscnoor [15] a montré que la teneur en gaz d'un 
tel systéme gaz-liquide augmente avec le débit de gaz 
jusqu’a un point critique, a partir duquel la teneur en 
gaz diminue pour augmenter de nouveau quand on 
continue 4 augmenter le débit de gaz. Une courbe don- 
nant la relation entre la vitesse ascendante des bulles de 
gaz et la teneur en gaz du systéme montre qu’on peut 
distinguer trois régions en augmentant le débit du 
gaz: d’abord, une région dans laquelle la teneur en 
gaz augmente tandis que la vitesse ascendante des 
bulles reste plus ou moins constante, ensuite une ré- 
gion ot la teneur en gaz reste plus ou moins constante 
tandis que la vitesse des bulles croit; finalement une 
région dans laquelle l'une et l'autre augmentent si- 
multanément. 


5. Travai sur Lit Mopite 
Dans les convertisseurs sur lit mobile on a un ensemble 
de particules ayant un contact direct entre elles et gar- 


dant approximativement une position fixe les unes par 


rapport aux autres, l'ensemble étant en mouvement 
par rapport aux parois du réacteur. Quoique pour ce 
mode opératoire on puisse s'attendre a rencontrer éga- 
lement des problémes hydrodynamiques et de trans- 
mission de chaleur, jusqu’A présent trés peu a été 
publié A ce sujet dans la littérature. En ce qui con- 
cerne la perte de charge dans un lit mobile nous 
pouvons nous référer 4 une communication de New- 
TON, Dunnam et Smupsown [12]. 


6. Travaw sur Lit FLvipist 
Un apercgu général des phénoménes caractéristiques 
d'un lit de catalyseur fluidisé a été donné par Ka- 
Bach [9] et par WiLHELM [18], qui considérent spé- 
cialement |’état de fluidisation dense. La plupart des 
études publiées traitent surtout de la perte de charge 
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dans le lit fluidisé dense: pour la fluidisation ag- 
grégative (WILHELM et Kwavk [19], Leva et collabo- 
rateurs [11]) et pour la fluidisation particulaire (W1L- 
HELM et Kwavk [19], Brtnkman [4], Verscnoor [16)}). 
C’est surtout pour cette derniére forme de fluidisation. 
qui se fait A l'aide d’un liquide, que le caractére de la 
relation entre la vitesse linéaire du liquide et la hauteur 
du lit dilaté, done la porosité correspondante, est 
établi. Si un lit de catalyseur est en état de fluidisa- 
tion particulaire , la porosité de la phase dense est 
uniforme. D’autre part, dans un lit en état de fluidi- 
sation aggrégative, on peut observer des bulles d’une 
grande porosité traversant un systéme d’une porosité 
inférieure. Le systéme présente alors l'image d’un 
liquide bouillant. Cette derniére forme est observée 
surtout si la fluidisation est faite A l'aide d’un gaz, et 
elle est caractéristique pour les lits fluidisés a l’échelle 
industrielle. A cause du brassage intense le contenu 
d'un réacteur sur lit fluidisé est presque homogéne, et 
par conséquent la transmission de chaleur du systéme 
aux parois du convertisseur est trés élevée de sorte 
que les gradients de température, qu'on signale souvent 
dans les convertisseurs sur lit fixe, ne se montrent 
pas dans les lits fluidisés. 

Au sujet des études de transfert de matiére, trés 
peu a été publié, tandis que pour la transmission de 
chaleur, VREEDENBERG [17] a donné un apercu de ses 
résultats ainsi que de ceux d'autres auteurs. VREEDEN- 
BERG trouve des coefficients de transmission de chaleur 
variant de 95 kealm?. °C. h. jusqu’é 450 keal m?. °C.h. 
selon les conditions de l’écoulement. 

Finalement il faut signaler la communication im- 
portante de Rets [13] qui, ayant rappelé les résultats 


des études pour le lit fluidisé dense, a présenté une 


étude sur le lit fluidisé léger, c’est A dire l'état d'une 


suspension dispersée. 


7. CONCLUSION 
Nous avons passé en revue un nombre important 
de problémes de Génie Chimique qui se posent dans 
la mise au point des procédés catalytiques. I] apparait 
qu’en ce qui concerne les problemes d’écoulement 
des fluides, de transfert de chaleur et de matiére, 
nous soyions en mesure de calculer un convertisseur 
a lit fixe. 
types de convertisseurs et il nous reste encore de 


Il n’en est pas de méme pour les autres 


nombreux problémes dont la solution demandera 
d’autres études approfondies. 
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NOTATION 
Ap _ perte de charge 
densité de fluide 
vélocité linéaire 
hauteur du lit 
surface extérieure du catalyseur par unité de volume 
du réacteur 
porosité du lit 
constante sans dimension 
eddy conductivity ou eddy diffusivity 
diamétre du catalyseur 
coefficient de conduction thermique 


chaleur spécifique 
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Summary, /'art / 
conductances to air flowing through a heated tube 


A preliminary study has been made of the effect of various end conditions on surface 


The effects of thermometric errors and flow disturbances lead to characteristic forms for the variation 
of surface conductance with flow charactermtics. The inlet disturbance produced large changes from the 
accepted correlation forms for nominally streamline flow for values of Grarrz number between 5 and 
300 in tubes with lengths 12 to 48 times their diameters. 

Part ll 
heated tube to air flowing through the tube of | in 
\ range of LD, of 4:1 has been tested. 
250 (L D,) &" 


The surface conductances for region (b) agre« 


Surface conductances have been measured in the transfer of heat from the wall of a steam 
diameter, the flow being disturbed at the inlet by a 
perforated plate and gauze Three regions of flow have been 
observed in the range of He from (a) 100 t (b> 20 D,) °** to 2200, fe) 2200 to 4000, 


with the equation 


L as 
lh 


(Re )0 70 


In the flow range wherem DG x 2200) the following equations apply 


Nv 0-05 (Gry? tor (fir) es 


tor ((ir) ys 


Résumé 7’ 


wedi ne nenat sp ren 


botuck 


ifique de transmission thermique d'une surface vis-a-vis d’air s’écoulant dans un tube chauffe 


emiere parti preliminane de Vinfluence de différentes conditions terminales sur le 


Les effets des erreurs thermore triquies et dles perturbations dans Pécoulement conduisent a des formes 
caractéristiques pour la variation de la conductance avec les caractéristiques de écoulement, La perturba- 
tron d’entree introduit des écarts considérables avec les correlations usuellement acceptéees pour des écoule 
ments lamimaires (nombre de Grarrz compris entre 
{8S cdiameétres 


Th urieme 


Jet 300, tubes dont les longwueurs varient entre 12 et 


partic. Mesures de conductances de surface dans léchange thermique entre la paroi ¢ hauffee 


& la vapeur et Pair traversant un tube de 25 mm diam. A lentrée, lécoulement est modifié par une plaque 


perforée et une toile. Les expeériences ont couvert une garmme de 4a Fig. | pour LD, et trois régimes d'écou- 


lement avee Re compris entre 
im et 250 (L/D, 
weil D,)®.*6 et 2200 


vor et 4000 


Pour le reritin b » les resultats sont been re pPreserites peer 


Vu OT6(L D, 


@,75 Ry 0)0 70 


tandis quavee 1,6 4 2200), les reactions deviennent 


Wr O45 0.70 pour itir) 


pour 


Part I. Effect of entry and exit conditions 


During work on the heat transfer from a heated 


The apparatus used is shown in Fig. 1. Steam 


tube to air flowing through the tube [4|, it was under a slight gauge-pressure was supplied to the 


The 


with the inlet and outlet temperatures of the air 


observed that slight modifications of the inlet and annular jacket. steam temperature together 


outlet design caused appreciable differences in the 


results obtained from the terminal 


temperature 
measurements. A variety of simple inlet and outlet 


devices were tested on a short heating tube, since 


the average conductances for long tubes are less 


sensitive to inlet conditions, in order to form an 


opinion on the importance of these points. 


stream were measured with thermometers. The air 
flow was measured with capillary flowmeters, which 
had been calibrated against standard meters and 
the results charted according to the method of 
Waurrwe.t {12}. The inner tube was fitted with a 


perforated plate covered with copper gauze of about 
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1} mesh to provide a support for packings such as C Perforated plate 
covered by wire Cc 
lead shot, glass beads ete., when testing the effect = gouze 
of such packing on the heat transfer from the tube Cc : 
wall to the gas. The heat transfer rates were es J 
timated from the flow rate and bulk temperature C ‘R R CG 
rise of the air stream passing through the tube t { 
oa : : . (a) (d) 
The surface conductance for heat transfer from the ¢ ¢ 
. or G = 
wall to the air was assumed to be equal to the over 
all coefficient from steam to air considering the C 
L = 956 to 317 in Cc 
Y/D, 7 64 fo 260 
ee 
f Steam —>= — rd 
\] Copy ' 
| 1 J Baffle N 
| N WN 
Nf 
NH 
| IN N 
NA 
7 Packed R N N LA 
3 siding NA 
Tus ‘ > 
| 2 Wa | ‘c) y 
Ame 
el | O-12 R 
| Plate and - te) 
| Gouze 
yu? | oo Rubber R 
Py ie -- Lo C 
s A evoder Cok Ch Inlet Conditions 
Glass CG Dp = 122 in 






Fig. 2. Inlet arrangements. 





Glass 





¥) Vent to atmosphere 






Fig. 1. Simplest heat transfer tube arrangement. 









steam condensate film and copper wall as offering 


negligible resistance to heat transfer compared with 








the resistance of the air film side. The area for 










ib) without cork insert 


heat transfer was taken to be the heated length 
ic) with cork insert 





of the steam jacket. The jacket consisted of a tele- 






scopic tube arrangement which was varied in length 





to suit the various inner-tube lengths. The heated 






length was varied from 31-7 to %56in. giving a 
LD, ratio from 26-0 to 7-84. 


woe ey 






The following discussion of the effects of inlet 






and outlet fittings is based mainly on tests of the 
L dD, 7-84 tube. 















Experimental Results (d) 
The effect of the inlet and outlet designs shown in Class C 
Fig. 2, 3, respectively were studied. The data ob- an ie f 
: , me s a or 
tained are shown in Fig. 4, 5, 6, and the key to Insulohion AL YT 
Rubber R 





the arrangements is given in Table 1. The experi- 






mental results have been based on the arithmetic Fig. 3. Outlet arrangements. 
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The effects of end conditions on heat 
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Fig. 5. The effects of end conditions on heat transfer coefficients. 


mean of the terminal temperature differences between 
the steam and air in order to provide ready compar- 
ison with the usual expression for the Grarrz func- 
tion and its associated equations. The GRAETZ equa- 
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transfer coefficients 


Table 1, Arrangement of inlet and outlet fittings 


for L Dd, 7-S4 





Perforated plate 
+ with 
without 


Inlet (see Outlet (see 


Run Fig. 2) Fig. 3) 


Figure 


~-— & & ee +e ee — 





tion for heat transfer to a stream having parabolic 
velocity distribution 


is given [9] as 
2,;We\/l 
NG 


h a l 7 
ile kL 


— = = (1) 


— eer) 
+ S@(n) 





Vol. I 
v - 1952 


where g(n) = 0-10238¢-M-8272" 4. ()-)1 220 > 89-22” 
0-00237 ¢- #2" . 

a (kL 
4 | We): 


This form is shown in Fig. 4, 5, 6, together with the 
approximations 


» 
N “ = (Gr) 


for Gr values up to about 3, and 


Nu = 1-75 (Gr)! 
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problems are not available. The above generalisation 
seems appropriate in view of the development of 
parabolic velocity distribution for laminar flow 
probably requiring [2] a transition length of 0-13 D,. Re 
The decay of turbulence may take a length equivalent 
This 


decay is also sensitive to the turbulence level at the 


to a few hundreds of grid-mesh widths [6]. 


inlet to the grid, which has been observed here and 
elsewhere [1] in the effect of inlet sections. 
The curves in Fig. 4, 5, 6 all agree with 


GRAETzZ function at low rates, Gr<3, where 





for Gr> 10, up to the limit of streamline flow 
in a particular tube. The tube of LD, = 7-84 
would give a Reynolds number of about 2000 
near Gr = 150, and the data at higher flow 
rates corresponding to the higher Gr values 
would be expected to show evidence of the 
beginning of transitional flow before attaining 
full turbulence at perhaps Re>10000. The 
disturbance of the inlet flow pattern caused 
by perforated plates, flow direction changes, 
expansions of cross-section to flow etc., is ex- 
pected to affect the surface conductance func- 
tions mainly in the transition range from the 
relatively undisturbed laminar motion at very 
low Gr up to the fully turbulent condition at 
high Re. 


The data have been grouped in Fig. 4, 5, 6 














CHOLLETTE__. 





Lu iil 4 


l 1 





using the data from Run 18 as the common 
curve as a guide to comparisons. 

The data for Run 32 in Fig.5 fall consider- 
ably lower in ordinate than the others, approach- 
ing the GRAETZ curve more nearly than any other 
arrangement. The approach to the GRAETZ curve is 
probably fortuitous rather than by the production of 
parabolic velocity distribution and thus conforming 
to the basic flow regime used in the GRAETzZ analysis. 
The data in Fig. 6 give curves similar to that observed 
by CHOLETTE [3]. 


Discussion of Results 


In all the experiments in Part I and Part II the 
heat transfer length was less than that probably 
required for the decay of eddies caused by the grid 
support and the establishment of the velocity di- 
stribution appropriate to the stable flow regime. 
Though studies have been made of the decay of 
turbulence induced by grids [6], [7] and of the 
development of parabolic velocity distribution from 
rod flow [2], [10], [11], solutions for the combined 


Fig. 6. 


~ 50 100 300 

kL 
The effects of end conditions on heat transfer coefficients. 
disturbances due to inlet fittings and errors due to 
radiation to thermometers have little effect and the 
heat transfer system gives a surface conductance 
relationship near to the limiting eq. (2). With Gr 
increasing the curve extends to higher ordinate 
than that on the Grarrz curve before the gradient 
falls to a value, which in many cases could be con- 
sidered as 0-4, agreeing with the CHOLETTE type of 
curve. Thereafter the curve changes to a gradient 
of 0-7 to 0-8 for the flow region up to and beyond 
Re = 2000, 

The experiments were not undertaken as a study 
of turbulence promoters. Two points were at issue; 
how much the surface conductance and the form of 
its variation with flow rate were affected by the 
differences between simple inlet devices, and how 
far the results were subject to thermometric radiation 
errors with various outlet designs. The major study 
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was to be made on packed tubes and thus the designs 
tested here were only considered qualitatively rather 
than with any intention of quantitative correlation 
or description of either the turbulence promotion or 
corrections. The form of the 


the thermometric 


hy 1M . 
kb? against the group 
~ were of interest for comparisons between the data 


4 / 
curves for NusSSELT number, 
We 


KL 
for packed tubes and for empty tubes with disturbed 


were 


outlet 


The variations in the outlet 
the 


inlet conditions. 


mainly to reduce radiation errors at 
thermometer, which could give estimated errors up 
to 12% in Nw for the outlet “a” on Fig. 3 and this 
error would be similar in packed and empty tubes 

The curve for Run 18, used as a “datum” for 
comparison, was provided by the simplest inlet (a) 
and outlet conditions (a). The combination of inlet 
(Run 25) had no effect on the 


lower 


“d” and outlet “ec” 
The inlets “b” and were similar and 
“a” the Fig. 4 


The omission of the perforated plate 


curve. 


with outlet gave curve on 


(Runs 22-27). 
and gauze cover moved the curve from the lower 
(Run 27) to the upper position (Run 28), when using 
the inlet ““d”’ and outlet “d’’. A similar but greater 
effect is noted on comparing the Run 30 (Fig. 4) 
with Runs 21, 31 (Fig. 5), wherein the former was 
without the perforated plate, the inlet “ec” being 
“d". The plate reduced the effect 
the baffle in the 


The glass evlinders used in inlets “b”™ 


used with outlet 
of the 
inlet “ec” 


disturbance initiated by 


“ce” and “d” were expected to reduce the heat 
transfer to the air from the metal wall below the 


perforated plate. The data for Runs 25, 


26 show the 
effect of including insulation in the head of the 
tube to reduce transfer to air from the wall outside 
the nominal heating length. The effect is small but 
not negligible. 

The use of the expanding inlet “e’’, seven dia 
meters in length, below the perforated plate had 
a pronounced effect, shown by comparing Runs 31, 32 
on Fig. 5. The omission of the perforated plate made 
Run 33 data differ considerably from the results of 
Run 32. As with Runs 30,21 the removal of the 
perforated plate lowered the flow rate at which the 
transfer coefficient curve passed to a gradient of 
about 0-75. The major difference between the curves 
for Runs 18, 33 was due to the reduction of the 
“jet” caused by the small tube in inlet “a”, and this 
effect was still further removed by the insertion of the 
perforated plate beyond the expanding inlet section. 
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The curve on Fig. 6 for Run 20 (inlet “ec”, out- 
let “*b”’) shows a difference from that for Runs 21, 31 


( . 


is mainly due to the greater radiation error in outlet 


on Fig. 5 which had an outlet This difference 
“d", despite the fact that both designs have high 
the 
radiation to the thermometers was best obtained by 


air velocities round thermometers. Reduced 


which insulated the thermometer 
the hot 


using the outlet “e” 
bulb 
provided high air velocity round the bulb, and also 


surroundings from metal surfaces, 


reduced the heat transfer from metal wall to air 
at any position higher than the limit of the nominal 
heated length at the head of the steam jacket. 
The inlet “e” and outlet “e’’ were adopted for 
further studies. It is not possible to discuss at length 
the effects 
Table 1, 


mind, e.g. diverting the air “jet” at the inlet, reduc- 


of each of the combinations listed in 
but each variant had a particular point in 


ing transfer from metal walls outside the nominal 
heating length, reducing radiation to thermometers 
measuring air temperatures. 

In comparing data from packed and empty tubes 
it will be necessary to assess the latter which are 
The 
opinion that heat transfer in packed tubes is greater 


likely to be used as a “datum” variation. 
than in the empty tube depends on the flow region 
and the degree of disturbance at the empty tube 
inlet. At Gr <3 it is very unlikely that the introduc. 
If the 
inlet conditions are similar to those in Run 18, the 


tion of packing has any appreciable effect. 


effect of packing may be small over a very wide 
flow range. 

Another source of error not mentioned above lies 
in the possibility of radiation to the inlet thermo- 
meter. The above arrangements do not suffer from 
this defect, the radiation to the thermometer stem 
being small and the air velocity around the bulb 
being high. That the error can be high, even when 
the radiating surfaces are at temperatures less than 
200° F, 
thermometer at the inlet so that the bulb is near 


is demonstrated readily by inserting the 
the perforated plate. Fig.7 includes data for Run 3a 
on a tube packed with lead shot with inlet “a” 
the 
thermometer was placed near the plate. The inlet 


and outlet “a”, and for Run3b when inlet 
temperature in Run 3b fell from 73 to 43°C as the 
flow rate G was increased from 14 to 540 Ib/sq ft hr. 
The inlet air temperature was 17 to 22°C in all 
the tests. This great error in the measurement of 


the air temperature is, however, less important than 
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first sight. It affects both the 
estimation of heat transfered from the air flow rate 


would appear at 
and temperature rise, and the measurement of the 
mean temperature difference between the heated 
wall and the air. The effect on the final plots is 
exemplified by the lines in Fig. 7 which also includes 
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The effect of thermometric error at the inlet 


of a packed bed. 


Fig. 7. 


a correlation equation proposed elsewhere [5] 
h ly 09/ Dy \ 


113) PG var 
| j | : 
Ty! " 


P L 
134 | D; | (4) 


The data in Fig.4 are based on the logarithmic 


mean of the terminal temperature differences bet ween 
the 


heated wall and the air stream. The eq. (4) 





The following experiments were undertaken to check 
the high power of mass flow rate associated with 
surface conductance found in earlier tests on a short 
tube with disturbed inlet [1] in the flow region 
below Re 
streamline. The effect of L D, ratio was tested using 


2000, wherein the flow is nominally 


the tube lengths prepared for the associated work 
on packed beds [5]. It was expected that the ratio 
would have an effect on the deviation of the heat 
transfer relation from the eq. (3) for streamline flow 
when Gr > 10. 


Part IL. Effeet of length to diameter ratio 
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was obtained with an inlet of type “‘e” and an out- 


let of type “e”. In Fig. 8 examples are given of 
data from Leva [8] for tubes packed with glass 
The data are for down- 
The 
stated to 


have been placed close to the packing and appears 


beads, and with lead shot. 


ward flow of air through the packing. inlet 


thermometer in Leva’s apparatus was 


to have had no radiation shield. The radiation error 


in Run 3b reduced the ordinates on Fig.7 by about 
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Fig. 8. Comparison with Leva’s data for packed beds. 


reason why Leva’s data in Fig. 8 fall below the 


correlation found in this laboratory 


The 


independent of packing conductivity. The Leva out- 


using more 


shielded thermometers. discrepancy is also 


let temperatures were also stated to have been 
measured by a thermometer close to the perforated 
plate support, and these measurements also may 


have been in error. 


Experimental Technique 
The apparatus used has been described [5], and used 
the inlet section “e” (Fig.2) and outlet section 
“e” (Fig. 3). 


was heated by condensing steam. The overall trans- 


A thin copper tube of Lin. diameter 


mission coefficient from steam to air was again assumed 
to be the surface conductance from the tube wall 
to the air flowing through the tube. The heat transfer 
rate was measured by the mass flow rate and terminal 
bulk temperatures of the air. The surface conduct- 
ance has been based on the arithmetic mean tem- 
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perature difference between the tube wall and the 
air. The physical properties of the air were assessed 
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Fig. 9. The dependence of surface conductance on flow 
rate for various tube lengths. 
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Fig. 10. The dependence of surface conductance on flow 
rate for various tube lengths. 


at the mean air temperature along the tube. The 
various lengths of heating surface were obtained by 
cutting the tube and its jacket, and re-making the 
tube outlet head [5). 
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Experimental Results 


Within the limited experimental range the data on 
Fig. 9, 10 show three regions wherein the transfer 
coefficient may be expressed as a power function of 
1hG 
it 
The region at Re > 2200 contains converging lines 


mass flow rate, or of 


for different L/D, ratios which appear to approach 
the gradient 0-8 commonly found in turbulent flow 
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Fig. 11. The analysis of the length effect. 


studies. In this intermediate region, preceding the 
establishment of fully turbulent motion at perhaps 
about Re — 10000, the converging lines suggest that 
the effect of the L, D, ratio may become less significant 
in the turbulent region. 


The lines for different L/D, ratios are parallel at 
Re < 2200. The break point between the region of 
gradient 0-7 and that of gradient near 1-0 is a func- 
tion of the ReyNo._ps number. From a logarithmic 
plot of the critical values, Re,, it was found that 
the limit may be expressed as 
L \ouw 
Di) 


Re, = 250 ( 





Voll 
No.5 


The region from Re 200 (LD,)*"*® up to 2200 


contains lines of gradient very near to 0-70, and 
having a pronounced L PD, function. The correlation 
equation for the data may be found by plotting 
Nu Re®? against LD, logarithmically as in Fig. 11. 


The data are correlated by the equation 


. = A 075 Fe 
Nu = 0-76! p, | (Re)? , (t) 
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gradient —0-9 from Fig. 11 the correlation becomes 


. = L 09 
Nu = 0-205 | Dp) (Re)! ®, 
t 


Since the data for low LD, ratios show a consider- 
ably lower power for ReyNoLps number the use of 
the eq. (7) may be too limited. In view of the approx- 


imate equality of powers for (1.D),) and Re a satis- 
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Fig. 12. 


This equation includes the tube diameter in the 
conventional arguments but its effect as an independ. 
ent variable has not been tested. It is to be expected 
that this equation will be limited in application to 
systems with the high degree of disturbance at the 
inlet similar to that induced by the present plate 
and gauze fitting. 

The data in the region of lowest flow rates show 
straight lines whose gradients increase with in- 
creasing LD, ratio, though this increase is more pro- 
nounced with the low ratios. For ratios between 
24 and 48 the gradient could be taken as 1-0 and an 


equation deduced similarly to eq. (6). Accepting the 


Comparison with correlations for fully developed flows. 


factory correlation is to be expected by incorporating 
the tube length in Gr. The data plotted as CC in 
Fig. 12 correlate satisfactorily on one line for all tube 
For 
Gr~> 5-8 the line has a gradient 0-70 corresponding 
250 (L/D,)%" 
-5-8 lie on a curve 


lengths over the whole experimental range. 


to the region in Fig. 9, 10 between Ry 
and Re — 2200. The data for Gr 
which tends with decreasing flow rate to become 
parallel to the eq. (2), though giving lower values 
of Nu at a specific Gr. The eq. (2) is the limiting 
form for low flow rates of the Grarrz function for 
a parabolic velocity distribution in the tube, given 


above as eq. (1) and shown in Fig. 12. 
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From the two straight lines CC on Fig. 12, cor 
The 


224) agree with the equation 


relation equations can be obtained. data for 


Gr>5S8 up to Re 


Nu OOS (Gye? | (S) 


The data for Gr<—5-S8 agree with 


Nu ODS (Gr). 


Discussion ol re sults 


The surface conductances measured were the average 


values for the tube surface from the entrance up 


to the stated LD, value 


Fig. 13 contains the data 
LD, 


together with an example from CHOLETTE’S results [3] 


plotted as 4, cG@ against Re for each ratio, 








13. The effect of length on mean surface conduct- 


ances 


from a tube of LD, 42. The latter values are 
those based on arithmetic mean temperature dif 
ference given in CHoLerte’s Table 4. The curve for 


LD, 


375, which corresponds to the divergence between 


42 crosses CHOLETTE’S curve at about Ry 


the two sets of data at Gr — 5-6. The two curves 
in Fig. 13 are similar in having a minimum at about 
Re SOOO. 

The variation in surface conductance with length 
of heating surface is shown in Fig. 14. These data 
may be compared with Cuoverre’s Fig. 4 for local 
coefficients A,. The coefficients 4, were measured 
for a short length of tube about the stated LD, 
position from the tube entrance. The average con- 
ductances reported here will continue to fall when 
h, has attained a constant value with increasing tube 


length. 
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On the plot of Nw against Gr on Fig.12, both 
the Cnuo_erre data and the present data show a 
rapid change of gradient around Gr = 5. The results 
above this value of abscissa correspond to a line 
whereas the CHoLETTE has a 


of gradient 0-7 line 


gradient of O4. A discrepancy has been noticed 





10 











L 

DY 

Fig. 14. The effect of length on mean surface conduct- 
ANnCces, 

between the curve CC and that obtained from an 

earlier but similar apparatus (Fig. 1) shown as DD 

1-22 in. 


L = %56in. Curve DD changes gradient abruptly 


on Fig. 12. This apparatus had D and 
25 to a value of 0-4 and then increases 
therefrom above Gr — 140-150. 
sponds to the transition region at Re >0200 for 
LD, =7-34. The Cuoverre data were all for Re < 2200, 
and hence do not show this increased gradient. The 


at about Gr 
This increase corre- 


present results for the curve CC appear to omit 
the 0-4 gradient region and show a transition region 
gradient of 0-7 for Gr>5-6, probably leading to a 
slightly higher gradient for the turbulent region. 
Had the data confirmed the CuoLerre or curve DD 
pattern they would have shown a change to gradient 
0-4 at Gr 
to a gradient near 0-7 at about Re 


5-6, followed by lines breaking away 
2000 for each 
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tube length. The present tubes of L/D, from 12 to 
about Gr from 96 to 24 
this 


seen in curve CC on Fig. 4, showing that the inlet 


48 would show breaks at 


respectively. No indication of form can be 
disturbance has eliminated the 0-4 gradient region 
and extended the higher gradient transition region 
from Gr — 5-6 upwards. That the curves will change 
gradient slightly from 0-7 near the upper limits for 
the GRarTz group mentioned above is more obvious 
from the plot in Fig. 9. 

The difference between curves CC, DD is marked, 
but not readily explained. The curve DD is for a 
shorter tube than those appropriate to the CC data, 
but the inlet and outlet conditions to the tube were 
the same in both series of tests. The only observed 
difference between the DD and CC tubes was that 
the DD tube 
made up with an expanding jacket and all joints 
The CC 


number of times during the fabrication of the unit 


was a smooth surface drawn tube 


with rubber bungs. tube was heated a 
and was cleaned with emery paper before use. The 
tube surface was considerably rougher than that of 
the DD tube, and this may have been the major cause 
of the discrepancy between the two curves in Fig. 12. 
The data CC for the tube D, lin. are similar 
in form to those given as 4A from the tube D — 1-22 in. 
using the inlet section “a” (Fig. 2) and outlet sec- 
tion “a” (Fig. 3). The outlet condition in the latter 
case was likely to cause thermometric error due to 
the radiation from the hot wall. The probable error 
has been calculated and gives the change from AA 
to data BB in Fig. 12. The change from BB to CC 
is thus probably due to the flow changes induced 
a” by the sections “e” 
The curve DD for the tube D, 1-22 in., L 


%6 in. is the only one tending towards the accepted 


by replacing the end sections © 


correlation for turbulent motion as given by the 


equation 
Ady 


0-023 ( Re y" | T 


(10) 


This equation is for surface conductances based on 
the logarithmic mean of the terminal temperature 
differences but the difference from /, should be 
The 
gradient 0-8 on Fig. 12 show the positions appropriate 


o 


within 5% for all the present data, lines of 
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to eq. (10) for various LP, ratios, the lines being 
extended down to Re — 2000, though they are not 
expected to be accurate at such low levels of turbu! 
ence. The equation approximates to the data to be 
expected for fully developed flow since it was obtained 
for long tubes. The data CC are for all tube lengths 
tested and indicate that the disturbed inlet conditions 
have little effect on the longer lengths but increase 
the surface conductances for the shorter lengths in 


the flow range tested. 


NOMENCLATURE 
specific heat of air at constant pressure; B.Th.U. lb F 


diameter of tube; ft (unless stated otherwise) 


diameter of packing; ft (unless stated otherwise) 


flow rate: lb/hr sq ft, (based on tube cross-sectional 


area) 


surface conductance; B.Th.U./sq ft hr °F; 4 based on 


logarithmic mean of terminal temperature dif 
ferences, h, based on arithmetic mean of terminal 
differences, h, surface con 


temperature local 


ductance 
thermal conductivity of air; B.Th.U.ft/sq ft hr °F 
length of heated surface; ft 
flow rate; lb/hr 


Viscosity of air; Ib ft hr 
We ‘ : 
. dimensionless 
KL 
}, dimensionless 


. dimensionless 
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Heat transfer to air flowing through packed tubes 
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(Received 6. December 1951) 


Summary —A study has been made of the transfer of heat from a heated wall to air flowing through packed 
beds in a tube of | in. diameter. The data are correlated by equations: 
, 1h G 
(a) tor law: a Hw 
Mu 
A Dy DD 113 L 0%”, Pig wi 
(134 ef} | 1—t 
k ly ly “ 
. Th 
(b) for 1600 S500; 
Mu 
A Dy : / Lae D, 0 sad DG 
5 
A 1 Ly 7 
IL 0165 
mm" then | 
ly 


Phe thermal conductivity of the packing has no influence in region (a) but affects the transfer in region (b). 
The flow rates used were lower than those reported by Leva ef alii ({ 1], [2], [3)]), but it appears likely that 
the present data trend to agreement with the functions proposed by Leva as the flow rate is increased. 


The most significant features of these results are that in the flow ranges studied the length of the packed 


tube is an important variable, but the packing diameter has little effeet on the heat transfer. 


Résumé — Etude expérimentale de la transmission de la chaleur entre la paroi chaude et lair qui parcourt 


un tube a garnissage de 25 mm de diameétre 


oxy rimentaunx 


Les equations ci-dessus (a et b) représentent les résultats 


La conduectibilité du garnissage est sans influence dans le régime a, mais affecte la transmission dans le 
regime b. Les debits utilisés sont plus faibles que coux de Leva et coll. [1], [2], [3]. Il) semble bien que les 


résultats actuels tendent vers les formules proposées par Leva, quand les débits croissent. 


Conclusion la plus significative a tirer des mesures 


dans la gamme des débits étudiés, la longueur du 


tube garniest une variable importante, tandis que le diamétre du garnissage n'a qu'un faible effet sur la 


chaleur transmose 


As part of the study of heat transfer in packed bed 
catalyst tubes, tests were made of the heat transfer 
from a heated tube wall to gas flowing through packing 
within the tube (Cau [4]). The primary object of 
these tests was to assess the effect of packing diameter 
to tube diameter ratio, packing length to tube dia- 
meter ratio and the effect of thermal conductivity of 
the packing material, as previous publications did not 
consider these factors (CoLBuRN [5], Leva [1 |). Owing 
to the limited apparatus available the flow rates were 
restricted to below 2000 Ib sqft hr based on the empty 
tube cross-sectional area. The packings available gave 
a 50:1 range of thermal conductivity. The L DP, ratio 
was varied over a 4:1 range. During the course of this 
work Leva ef alii ({2}, [3}) published data from sim- 
ilar studies. Their ranges of flow rate were mainly 
higher than those reported below, and it is considered 
that the present data are complementary to those of 
Leva. This consideration has been included with the 


discussion of the results. 


EXPERIMENTAL TECHNIQUE 
The apparatus used is shown in Fig.1 incorporating 


the inlet and outlet sections found most suitable after 


trials of various designs to minimise the errors in mea- 
surement of terminal temperatures for the air stream. 
The coned inlet section reduced the distortion of flow 
distribution to the packing and provided insulation to 
reduce heat transfer to the air from the inlet tube 
below the flange. The cork insulation at the outlet was 
drilled as shown to exclude radiation from the packed 
bed to the thermometer measuring the gas tempera- 
ture. The Lin. bore copper tube and the 245 in. bore 
jacket were reduced in stages from the length shown 
in Fig. 1 to give heated lengths of 48, 42, 36, 30, 24, 
I2 in. After each reduction in length the connections 
to the inner tube and the dimensions of the outlet head 
were re-made to those in Fig. 1. In each case the 
copper tube was filled with packing to the level of the 
top of the steam jacket, the effective length L quoted 
below being the steam-heated length. The cork 
inserts at the inlet and outlet reduced the error pro- 
duced by heat transfer from the tube walls outside 
the nominal heating length. The packings used were 
all either spherical or nearly so, and of smooth surface. 
The Socony-Vacuum (“Sova”) catalyst beads were 
not uniform in size. The mean diameter of a large 
number was 0-145 in. The glass beads were sufficiently 
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near uniformity to consider them of constant dia- 
0-191 and 
0-0388 in. The lead shot were uniform in diameter and 
used of D, = 0-255, 0-182, 0-126, 
00-0951, 0-0426 in. Polished 
steel balls of diameter 
0-250 in. 


The packing was sup- 
I } 


meter and packings were used with D), 


packings were 


were also used. 


ported on a perforated plate 
drilled with 4 in. dia, holes 
spaced fy in. apart. A dise 
of 60 mesh gauze was placed 
between the packing and 
the perforated plate. 

The air to the heated 


tube was supplied through 





capillary flowmeters cali- 
the 
technique of Wiuirwe  [6}. 
The 


sufficiently to reduce press- 


brated by plotting 


flow was smoothed 
ure fluctuations at the inlet 
to the bed to within 1 mm 
in 700 mm _ water 


cvauve 
ua uge 


pressure difference across 
the packed bed. The tempe- 
rature and pressure of the 
air were measured at stages 
the 


supply tubes for estimation 


through meters and 








o=4p of the mass flow rate and 


P =Manometer connection 
C = Condensete 
Fig. 1. Heat transfer 
apparatus, 


frictional energy loss in the 
bed. Steam temperatures 
were measured with thermo- 
meters at the inlet and con- 
densate outlet to the steam jacket, since sufficient 
excess steam was used to blow through the vents to 
remove rapidly noncondensible gases and condensate, 
and the pressure in the steam space was greater than 
atmospheric. 

The temperatures of the air stream at inlet and 
outlet to the packed bed were measured with thermo- 
meters graduated to 0-1° C and appropriately calib- 
rated. Even in this small apparatus tests had to be 
made for times up to an hour at the low air rates for 
ateady state conditions to be attained. 


EXPERIMENTAL RESULTS 


Tests were made for various air rates covering the range 
available from the existing equipment in the laborat- 
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Fig. 2. The effect of packing on the relationship betwee 
h De DG 
and Reynolds number 
k } , 
30 in. L 
1 in. D, 


Nusselt number 


Packing length 
tube diameter 


Packing Symbol Run 


lead shot 


dD, in. 
W255 
182 
126 
OOds 
O426 
191 
(-O388 
250 
145 


vlass S| sheres 


steel balls 


sova beads 
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Ly, Symbol Run Symbol Run 
48 @ 43 : A 54 
42 Oo 47 © 64 
36 A 51 : © 6S 


Glass spheres Dy, 0-91 in. 


ory. The series of packings were used for tests of dif- 
ferent tube lengths to assess the effect of the ratios of 
tube diameter to packing diameter and of the tube 
length to tube diameter. The results of these tests are 
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exemplified in Fig.2, 3, 4, 5. The data in Fig.2 show 
the effect of various packing diameters using a fixed 
tube length, while those in Fig. 3, 4, 5 are for one pack- 
ing and various lengths of bed. The surface conduct- 
ances reported for transfer in packed beds were based 
on the mass flow rate and mean temperature changes 
of the air stream, the tube inner surface, and the 
logarithmic mean of the terminal temperature differ- 


ences between the steam and air streams. The thermal 








tio 
ou 


resistances for the steam side and the copper tube 
wall were negligible in comparison with that of the 
air stream and the temperature of the inner face of 
the copper tube could be taken as that of the steam 
without appreciable error. The values of the physical 
properties of the gas used in graphs and correlations 
were assessed for the mean temperature of the air in 
the packed bed. 

The data all show lines of the same gradient for 


flow conditions below a specific Reynolds Number, 
DG , a 4s 
r . for each packing and tube length. The critical 
“a 
D@ ., 
value of : is 
“a 
(Fig.2) but is almost independent of packed length 
(Fig. 3, 4, 5). 


limit dividing two possible correlation regions of the 


affected by packing dimension 


As this breakpoint was to be used as a 


present data, it was of value either to define its 
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variation or to express it in constant terms. It was 
found that the variation in the critical Reynolds 
D,G@ 
it 
Fig. 3, 340-460 in Fig. 4, 390-450 in Fig.5 was reduced 


considerably by using the analogous expression for 


Number, of 56-400 in Fig. 2, 285-350 in 


the tube diameter. The groups above gave mean 
Db fi *: . ms 

values for ~"" of 1660 for the Fig.2, 3, 1560 for 

Fig.4 and 1610 for Fig.5, most of the values lying 


* 











8100 2 
fu —~ 
Fig. 5. 
Run D, Dn, LY 
e 440 255 ‘ 66 
40 O-255 l: 70 


Nymbol Symbol Run dD, dD, 
W255 
O- 255 
mOOD 
(L104 


(h255 Leva 


Le a i shot 


Fig. 3,4 and 5. The effect of length of heat transfer surface 


on the variation of surface conductance with flow rate. 


A value of 1600 was 
accepted for all packings and tube lengths as a critical 


well within + 10% of the mean. 


value of being the upper limit of the correlation 


region at low flow rates. 


Correlation for 1600 
The data in Fig.2, 3, 4, 5 show a common relationship 


A Ds 


D, G Wiz 
foe 
k mn" Mu 


(1) 


where « is probably a function of D, D, and of LD. 
The data for L — 30 in. 
function 


in Fig.6 agree with the 


(2) 
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and the data for the other lengths have similar fune- 
tions with 8 as a probable function of L/D,. 

The value § is shown in Fig.6a, to conform to 

0-90 

(3) 

The value of y is independent of the other packing 

characteristics, particularly thermal conductivity, the 

data for all packings and tube lengths tested for 

ynG 


1600 being in close agreement with the equation 


113, J 
D, | 


h 1h 
ik 


09, P,GAiw 
——"F (4) 


“u 


ose (7? | 
t 


which has thus correlated data for the ranges; D,, D, 
from 0-0388 to 0-255 and LD, from 12 to 48 based on 


BL 


Packing 
@ lead shot 


glass spheres 


Symbol 


A ateel balls 
© sova beads 














on - LAG 
big. 6. The effect of packing diameter when ae 
Mu 


1600. 
the | in. bore tube. It will be noted that in the deriva- 
tion of eq. (4), as with eq.(S), the dimension D, has 
been introduced in the conventional ratios D,, D, 
and LD, but has not been tested as an independent 
The work of Leva et alii ({1}, [2], [3}) 
suggests that correlations of the forms used here will 
cover the effect of D,. 


variable. 


hG 


Correlation for 1600 < 3500 


Unfortunately the equipment available limited the 
: ° hG - 
air tlow rate to that corresponding to ~‘ 3500, 
but the range of data appears to be sufficient to 
justify the following correlation. 


hG 
For 


line in Fig. 2, 3,4, 5 the gradient of the plots decreases 


> 1600, i.e. above the break point on each 


below 1-17. On Fig. 2 it seems that this lower gradient 
remains almost constant for various particle sizes at 
constant LD, ratio, though on Fig.3, 4, 5 the gradients 
vary considerably for different L/D, ratios. With the 
limited data available (Table 1) the gradient may be 
assumed dependent only on the length of the packed 
beds. The plot in Fig.7 shows the straight line drawn 
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through the values for glass beads, and this line is in 
reasonable agreement with the data for the other 





0-02 
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hG 


Fig. Ga. The effect of heating length when 1600. 


materials. The lines on Fig.2, 3, 4, 5, thus suggest the 


form 


h Ih D,G m 


at 
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The effect of heating length on the function of 
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S500, 


Table 1. Values of gradient m. Data from Fig. 3, 4, 4 





Lead shot 


0-255 in. 


Glass beads Steel balls 


0-191 in. 0-250 in. D 


LD, , 
Mv M P 


1-00 1-03 
1-03 
O-S9 
1-09 
1-07 
0-82 


1-05 
1-02 
0-96 
0-93 
0-90 
O-SI 


1-00 
0-96 
O-s82 





Data from Fig. 2 
All packing materials. LD, 





D, Dy, Gradient m 
1-05 
1-00 
0-94 
1-02 
1-03 
1-02 
O-95 
0-99 
1-02 


0-255 
0-250 
0-191 
0-182 
0-145 
0-126 
0-095 
0-0426 
0-0388 
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From Fig. 7, 


ie L ows 
ODD | De) (6) 


The data from Fig.2 have been analysed to find the 
value A, in eq.(5) using the power m from eq. (6). 


Pa hin ] ‘ rbd 


yvlass balls 

steel balls a 
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Sova beads 


oor 4 
0 or 
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Fig. 8. The effeet of packing dimension when 1600 


nwa 
a S500, 
fé 


The plot in Fig.8 shows a straight line for LD, — 30 


of gradient 0-90. Similar lines would probably 
describe the variations of A, with D, D), for other 
LD, ratios, 


i 0-20 





yo 6S 


#/ 





tno 
(2%, 





‘with m=0-$5 


| 
| 
mao 


“p “/u) 


9 Z> 
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20 2 ww 
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[(AOr/) (Op b, 





60 8 100 


hG 


» effect of heating length when 1600 
SOO). 


Fig. 9. T 
“ 


Thus accepting 


(7) 
is shown that the effect of thermal conductivity is 


K, is probably a function of L/D,. The data shown 
in Fig.9 are the mean values of A, for the values at 
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each LD, ratio. They conform to a straight line of 


gradient 1-82. The data for air flowing through 
various materials, packings and packed lengths are 


LhG 


all correlated for 1600. 3500 by the eq. (8), (Sa): 


Aly 182 Dy 0 90 DG m 
i : Dy! “ 
onl & 


m O90 | Dy } 


(8) 


0 165 


(Sa) 


. * ,c, 
A function of Prandtl Number | ) has not been 


included in the correlation equations. In view of the 
probability that the appropriate function varies less 
»o 


than 2% for the experimental range, such an inclusion 


cannot be assessed from the present data. 


Discussion oF ResuLTs 
It had been expected that the thermal conductivity 


of the packing would influence the heat transfer, and 


*» > 


* 
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Packing Symbol 

191 
-O388 
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W255 glass balls 


(Is2 
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0426 


The effect of flow rate on surface conductance for 
various packings and heating lengths. 


during this work a publication by Leva and Grum- 
MER [2] discussed results showing that the transfer 
coefficients for a similar system increased with thermal 
In Fig. 10 it 
ha 

. 1600, 


(} = 900 Ib hrsq ft, the 
effect seems to become apparent at higher flow rates. 
The data of Leva and Grummer [2] were obtained 


conductivity of the packing material. 
insignificant below the break point at 


approximately at though 
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at higher flow rates except for their tests on lead shot 
and aloxite particles. These tests are referred to later 
as the present data explain what Leva considered to 
be anomalous results for lead shot. Considering the 
48 in. and L = 42in. in Fig. 10, the 
lead Dd, 0-255 in., steel balls 
O-191 in. all lie 
For other 


lengths L 
points for shot 
D,, — 0-250 in., and glass beads D,, 
on one line for a given packed length. 
lengths of packed beds the data for packing of greatly 
differing thermal conductivity also lie near the same 
line. The effect of D, differences should be taken into 
account in this comparison but from eq. (4) it is seen 
that the effect of D, on h is small since hx D°", 
Thus a ten-fold change in PD, alters h by 10%, and 
the maximum ratio of two particle sizes in the present 
tests was only 6-6. The small effect of particle size 
in this range is demonstrated by the data in Fig. 10. 


LhG 


Above the limit near @ ), 1600. the 


7 
extent of the present data is sufficient to suggest that 
thermal conductivity of the packing begins to affect h. 
In Fig.10 the data for L = 24 in. and L 


show that at the highest flow rates, the lines for glass 


beads are falling away from those for lead shot and 
steel balls. This trend is also obvious from the gra- 
dients shown in Table 1. The present data thus appear 
to lead towards the results of Leva and Grummer 
showing the effect of thermal conductivity in the 
higher flow ranges of G — 500 to 4000, 

Leva and Grummer [2] found / proportional to 


DG 
i 


the OS to 9 power of except for their data 


on lead shot which have been included in Fig. 5. 
The 


to the lower flow rate region in the present work, 


Leva results were not anomalous, conforming 
probably being below the break point corresponding 
to those shown on Fig.5, whilst the rest of Leva’s 
results were obtained for the higher flow region 
where the power 0-S-0-9 would be expected from the 
present data. The existence of the region of high 
gradient 1-17 does not appear to have been recognised 
elsewhere in the case of fixed beds except by Scnvu- 
MACHER [7]. The present data confirm work on the 
fine powder packings (AGarwaL and Srorrow [8}) 
whether static or fluidized. The Leva [1] data for 
0-172, D, the 
trend to a higher power than 0-8 with decreasing flow 
The data of Leva, Wernrraves 
MER [| for fluidized beds appear to confirm the high 
gradient above. hx GY 


it has been shown (AGARWAL and Strorrow [8]) that 


glass beads D,, 2-0in. also show 


and GruM- 


rate. 


Their results show and 
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the appropriate plots of data for fluidized beds 
continue the curves for the associated packed bed 
at lower flow rates, and show the same dependence 
on mass velocity. The Leva power 1-15 is thus a pos- 
sible confirmation of the present power 1-17. There isa 
difference which must be noted in that the Leva, 
Werxtraus and GRUMMER surface conductances were 
obtained from the integrated mean temperature dif- 
ference based on a temperature survey along the flow 
path of the gas, whereas the present data are based 
on the logarithmic mean of the terminal temperature 
differences. 

The equations obtained above have been checked 
against the data used in their construction and it was 
found that the eq.(4) agreed with all experimental 
results to less than | 10% and with the majority of 
results to less than 45%. A comparison has also 
been made with the few data from Leva’'s tests which 
(a) those for glass 
beads 1D, (-172in. in a tube of D, 1-96 and 
L 36in., and D, the tube, 
from Leva [1], (b) those for lead shot of dD, 
in a tube of L I4in. and D), = 0-824 in. 


cases the eq.(4) gave A values 19-9 to 22-3% higher 


are in the same flow region. i. 


0-388 in. in same 
O-O9L in, 
In these 
than the observed values. On checking data for an 


earlier apparatus for thermometer arrangements 
similar to those employed by Leva it was found that 
they were about 14% higher than the eq.(4) values. 
The data referred to were considered unsatisfactory 
primarily owing to the radiation errors to the thermo- 
meters concerned, The data of Leva [1] and Leva 
and Grummer [2] may be in error by a similar 
amount. This may have contributed to the fact that 
on reversing the direction of heat transfer in their ap- 
paratus Leva, Weinrraus, Grummer and CLARK [3 | 
found surface conductances for cooling the gas about 
15% higher than those from the Leva correlations 
of conductances when heating the gas. 

Another difference between the present results 
and those of Leva lies in the function of the ratio 
D, D, adopted. Leva {1} used the exponential func- 


tion « ° Pp! 


’* to cover the wall effect dependent on 
D, D,. The data for the higher ranges of flow tested 
by Leva, for which the appropriate power of Rey- 
nolds Number was 0-0, agreed with a straight line 


when plotted as 
DG O98 


Aly 
k 7 


k "uy 


against D, D, as in Fig.5 of Leva [1}. The data for 
a similar range of D,, D, to that for the present results 
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could be approximated by a funetion 


hi 5; Dy G oe D, 
~ 


i ” lh 


similar to that used here, though the agreement be- 


tween the experimental points and this form is less 





Packing 
lead shot @ 


glass balls 
steels ball« a 
sova beads © 


‘ gue , 








0095 


‘ > 49 > 


EFT I 


Fig. 11. The failure of the exponential function for the 


effect of D, D,. 
than with the exponential form above. The exponen 
tial form may be correct for the higher flow rates but 


for the two regions in the present work, the form is 


A 
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Fig. 12. 


unsatisfactory for correlation purposes as shown in 
Fig. 11. Using the ordinate appropriate to each flow 
region the curvature of the plots is pronounced, 
though it should be noted that the data in Fig.5 of 
Leva [1] show similar curvature, though less pro- 
nounced. The simple power functions of D,, ), used in 
eq. (4), (8) will not show a maximum in a plot of 
The 


exponential function of D, D, used by Leva does 


AG” or h@" respectively against D, D,. 


show this maximum at J), D, — 0-15 in agreement 
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packing material, 
Various L D,, 
Constant D, dD, 
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with the experimental data at the higher flow rates. 
Leva, Weintravs, GrumMER and CLARK [3] discuss 
the use of these 2), D, functions. 


CONCLUSIONS 

From consideration of the present data together with 
those of Leva it appears likely that the heat transfer 
to gas flowing through packed beds may be covered 
by the three types of curves shown in Fig. 12. These 
types cover the probable effects of length of packed 
bed, the diameter of packing particle and the thermal 
conductivity of the packing material. 

The curves expressing the heat transfer to air 
flowing through packed tubes having the same ratio 
DL, and different ratios D, D, are of the form in 
Fig. 12a. The regions ab and bed are exemplified by 
the data given in Fig.2 for five sizes of lead shot for 
DL, = 30. the 
relative positions of the lines but do not alter the 


The ordinates used in Fig.2 move 
xz 


forms from those to be expected on plotting A instead 
of =. The lines for different D, D, ratios in 
Fig. l2a are parallel to each other in each of the 
regions ah, bed. 

The curves for heat transfer when the J), ), ratio 


is constant and LD, is varied are exemplified by 





y? 
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Different packing material. 
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The suggested ranges of effect of L, D,/D,, and thermal conductivity of the packing. 


- — . — ™ DG 
Fig. 12b, similar to the data in Fig.3, 4.5. At ~‘ 


“a 


above the critical point, the lines in be vary in gradient. 


This gradient increases with increasing L D, ratio. 
It is reasonable to expect that the lines will converge 
to give a common line ed. In this flow range ed as 
covered by Leva’s data the surface conductance is 
independent of LD,. Leva, Grummer, Wery- 
TRAUB and CLARK tested ratios LD, — 10 and 17 with 


DG d : 
: 640-1400 and found no length effect in the 
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cooling of gases flowing through packed beds. These 
data may not have extended down to the region ed 
in Fig.12b but the present data show that such 
LD, changes under these flow rates would have small 
effect on the heat transfer equipment. 

The vertical broken line in Fig. 12b cutting the 


lines indicates the critical flow condition (G — 9001b 


sqfthr: MO 1600 for D, — lin.) which depends 


on the ratio D, D, but not on L D,. 
the 
packed tubes having the same LD, and D,, D, ratios 


In Fig. l2e surface conductances are for 
but packings of different thermal conductivities. The 
thermal conductivity has negligible influence in the 
region ab but affects the heat transfer in the region 
he, leading to the effects measured by Leva and 
GrumMer [2] at higher flow rates, ed. 

The small effect of D, leads to a reasonable cor- 
Aly 
k 
observations of 


relation of all data on a logarithmic plot of 


1nhG . , 
: agreeing with the 
(6). This 


packings of various thermal conductivities and dia- 


against 


SCHUMACHER author commented that 


meters had little effect on the heat transfer between 
the retaining wall and the gas flowing through the 


packing. He found that below a critical Nusselt 


Aly oor 


-=), 


number, this group was proportional to 
DG a ori AD D, G 0% 

: _ and above this critical value“) '«({""’ "|. 
“ / k “a 

h Dy AD: 


57, the form ik 


Below a critical Nusselt number 


D, Ga ' . 
— } was found. SCHUMACHER considers that the 


~| 
Leva correlations based on | 
h ly 
k 
low flow ranges equivalent to the two lower regions 


D,G 0% rr 
*) are only valid in 


the region 57 - 225. The present data are for the 


of Scuumacner. The data in Fig.2 show an apparent 


— h 1% © -—— ry 
critical value of at about 30, rather than 57. The 
k 
gradients below and above their critical value agree 
with ScHuMACHER’S opinion. The lower critical Nusselt 


number corresponds to the lower critical value of 


DG ; 
" =: 1600 found here as compared to the value 


2300 found by SCHUMACHER. SCHUMACHER comments 
on the necessity for altering Reynolds number by a 
factor of 2 in the data of Kime [10] to make them 
agree with those of Leva, mentioning the 3:1 length 
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ratio though D), and D, were the same. The present 
analysis of the length effect explains the origin of 
such differences. 


Note added in proof— Papers in which heat transfer to fluids in 
packed beds has been analysed in terms of effective thermal con- 
ductivity and in which attempts have been made to segregate 
the thermal resistances in the system are intentionally omit- 
ted. These analyses, as exemplified by References [11 ]—[ 18], 
are based on various suppositions which are not correct. The 
present report is limited to a presentation of data in the form 
of surface conductances correlated by the usual arguments, 
rather than discussing the postulated mechanism of the heat 
transfer, The functions shown prove simpler than one might 
imagine in view of the probable complexity of the transfer 
system. It is not expected that the correlation equations for 
Re 1600 can be extrapolated below the experimental range. 
The values quoted are already lower than the constant value 
to be expected at low flow rates. The suggestion that the 
failure to attain the constant range at higher / values than 
those shown is due to heat losses from the apparatus [16 
The 


surface conductance and derived Nusselt number depend on 


seems questionable. numerical value of the limiting 
the basic definitions of the quantities involved, notably tem- 


perature distribution and effective conductivity. 


NOTATION 
specific heat of air at constant pressure; B.Th.U/Ib°F. 
diameter; D), for tube; ), for packing; ft. 
flow rate based on empty tube; lb sqfthr. 
surface conductance; B.Th.U/sqfthr’ F. 
thermal conductivity of air; B.Th.U ft/sqfthr F. 
length of heated surface; ft. 


viscosity of air; lb ft hr. 


REFERENCES 

(1) Leva, M.; Ind. Eng. Chem. 1947 39 857. [2] Leva, M. 
and Grummer, M.; Ind. Eng. Chem. 1948 40415. [3] Le- 
va, M., Weinrraus, M., Grummer, M. and Crark, E. L.; 
Ind. Eng. Chem. 1948 40747. [4] Cau, Y.C.; M.Se. Tech. 
Thesis 1949 (Manchester). [5] Cotpurn, A. P.; Ind. Eng. 
Chem. 1931 23 910. [6] Warrweti, J.C.; Ind. Eng. Chem. 
1938 30 1157. [7] Scuumacuer, R.; Erdél und Kohle 1949 
5 189. [8] Acarwart, O.P., and Srorrow, J. A.; 
Ind. 1951 (April) 321. [9] Leva, M., Wervrravus, M., Grom- 
mer, M.; Chem. Eng. Prog. 1949 45 563. [10] Kutne, G.; 
V.D.1.-Forschung 1938982. [ll] Brorz, W.; Chem. Ing. 
Technik 1951 23 408. [12] Copertey, C. A. and MarsHa.s, 
W. R.; Chem. Eng. Prog. 1951 47 141. [13] Hover, J. 0.; 
and Prret, E. L.; Chem. Eng. Prog. 1951 47 295. [14] Smrru, 
J.M., Moraes, M. and Spiny, C.W.; Ind. Eng. Chem. 
1951 48 225. [15] Smrre, J.M., Bunnews, D.G.; Irvin, 
H. B. and Otsox, R.W.; Ind. 1949 41 1977. 
[16] Verscnoor, H. and Scnurt, G.C. A.; Appl. Sci. Res. 
1950 A297. [17] Wiewe tm, R.H. and Sincer, E.; Chem. 
Eng. Prog. 1950 46 343. [18] Witwe om, R. H. and Bernarp, 
R. A.; 1950 46 233. 


Chem. 


Eng. Chem. 


Chem. Eng. Prog. 





D. W. van Kreveven and H. A. G. Cuermin: Erratum — Estimation of the free enthalpy 


Erratum: Estimation of the free enthalpy (Gibbs free energy) of formation of organic 
compounds from group contributions*) 


D. W. vAN KReVELEN and H. A. G. CuerRMIN 


Staatsmijnen in Limburg, Central Laboratory, Geleen, The Netherlands 


In this paper a number of numerical errors appeared. In pounds. The corrected values are given in Tables 13, 14, 16 
calculations involving thermodynamic properties of compounds — and 17 below (numbering as in the above mentioned paper). 
containing sulphur, bromine and iodine one is foreed to choose The Figs. 4, 6 and 7 should be changed in accordance with 
the gaseous state as the reference state. The .1//?-respectively — the corrected values, 


14?-values used in these calculations (see section V of our 


Table 14 





paper) should consequently refer to the gaseous state of 


B 
Group contribution A 4 T Highest 
, 1 
done in our paper. All values of the free enthalpy of formation tem perature 


sulphur, bromine and iodine. This has, erroneously, not been 


(iroup 


should for that reason be corrected by the free enthalpy of 300. 600K 600 > 1500°K of lit, value 


transition due to the transition of the elements from the Lewis ' B B 
, , 
and RANDALL reference state to the gaseous reference state. 


We recalculated the group contributions and free enthalpies 
(254 13-10 158 


0-26 1-62 0-26 


* Chemical Engineering Science 1051 1 66. 0-225 1-718 0-176 


of formation of sulphur, bromine and iodine containing com 


0 - 780 0 
Table 13 








B a” 
Group contribution A ae 7 Highest Table 17 





te mh pe reatlure 


0° 6600 OK 600° 1500 K of lit. value , 7" Cal. 
° ° Free enthalpy of formation A 
: culated 


B B : from 

00 * 600 K 600° aK values 

K qeven 

20-552 14026 21-366 «1-167 B by 
10-68 107 10-68 17 

3-32 1-42 3-32 1-44 SO, 86-657 +1-735 86-657 +1-735 L500 16 

OT = 0-51 065 O44 SO, 100-313 43-43 100-813 3-43 1200 51 

30-19 3-39 30-19 339 49-414 0-245 49-269 0-265 1500 1} 

82-58 55s sO-oo = 5-26 ‘Ss, 3-30] 0-150 3-104 0-177 1500 il 





Table 16 





0 -*#K 600 K 900 K 1200 K 1500 K 
Cal- 


Cal . Cal . Cal. . c Cal. . 
Compound . Found Found 7 Found al Found al Found culated 
culated culated culated culated culated from 


keal keal keal keal keal keal keal keal keal keal ‘Literature 


mol mol mol mol mol mol mol mol mol mol 


Methylbromide . ; 6-0 631 4+05 +00 4+ 76 + 7: +14 + 14-2 
Methyltetrabromide —. 0 + 69 - 16-6 
Methylbromoacet vlene 444 4443 “f 48-0 +514 
Vinylbromide . . . . +163 +154 4 +19-9 +244 
Acetylene di-iodide. . +669 466-7 4 +647 +610 
Vinyliodide . . . . . +265 4266 431-2 +31: + 36-1 
Ethylmercaptane Q- 9-5 

Dimethylsulphide . . > 70 

Thiophene .... . + 20- +20-6 +204 + 29-9 
Dimethylsulphoxide . 28: 28-6 ie 49 +196 
Dimethylsulphone . . 74-0 3: 43-6 12-9 
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K. G. Densicu; The Thermodynamics of the Steady State. 
Methuen & Co. Ltd., London 1951. VIII 103 pp. 6s. 6d. 


The Second Law of Thermodynamics is one of the most 
powerful tools of science when applied to equilibrium pheno- 
mena and reversible processes; in regard to irreversible pro- 
cesses, however, it can provide us with no more than the 
qualitative information that the entropy increases. The first 
attempt to extend thermodynamic methods beyond the 
equilibrium state was made by Lord Kevin not long after 
the formulation of the Second Law itself. Ketvin considered 
a thermo-electric circuit in its steady state, so that a constant 
temperature difference was maintained between the junctions, 
and a steady electric current flowed in the circuit. He assumed 
that, although irreversible processes were an essential feature 
of the system, the Second Law could nevertheless be applied 
to the parts of the system which were reversible, and he was 
thus able to derive relations between the magnitude of the 
Pevrier and THomson heats, and the temperature coefficient 
of the thermo-electric E.M.F. 
his treatment involved a new hypothesis in addition to the 


KELVIN realized clearly that 


Second Law, and stressed the need for experimental veri- 
fication of the resulting conclusions; subsequent data were 
certainly in reasonable agreement with the new equations. 

Later writers applied the THomMson hypothesis to a number 
of other problems such as the E.M.F. of a concentration cell, 
and the process of thermal diffusion (the change in concentra- 
tion associated with a temperature gradient). Irreversible 
processes were again characteristic of the systems, and the 
theory was applied to the steady state. For some of these 
problems a treatment based on the kinetic theory of gases 
was already available, and the results obtained by the two 
methods were in agreement. 

The major step towards clarifying the nature of the THom- 
8ON hypothesis was made by Onsacer in 1931. Chemists 
had always assumed that in an equilibrium involving several 
independent chemical reactions, there was detailed balancing 
of the individual reactions, and this was usually taken to be 
a consequence of the Second Law, Onsacer showed that 
although this was correct for a system involving two indepen- 
dent reactions, there was nothing in the Second Law to prevent 
equilibrium being maintained by cyclic processes when there 
were more than two independent reactions. A new hypothesis 
was required if this possibility were to be excluded, the prin- 
ciple of microscopic reversibility, and ONSAGER was able to 
show that THomson’s hypothesis was a consequence of this 
principle. ToL~MAN has expressed this principle as follows : 
“Under equilibrium conditions any molecular process and 
the reverse of that process will be taking place on the average 
at the same rate.” 

The past few years have witnessed a considerable develop- 
ment of the thermodynamics of irreversible processes, and 
a wide literature has sprung into being. Much of this litera- 
ture has been somewhat abstruse and scattered, and the 
subject has remained in the hands of specialists, and has not 
penetrated to the ordinary research worker. Dr, Dennicn 
has performed an extremely valuable service by making 
available to the wider scientific public a clear and readable 
introduction to steady state thermodynamics, and this mono- 
graph should go a long way towards dispelling the obscurity 
which has been associated with the subject. The treatment is 
elementary, and anyone with a knowledge of classical thermo- 
dynamics should have little difficulty in following the argu- 
ments. 

The monograph consists of 6 chapters. Chapter I contains 


a brief survey of the historical development. Chapter IT dis- 


cusses the application of the THomson hypothesis to thermal- 
migration phenomena, Chapter IIT is devoted to a discussion 
of the Onsacer theory. ONSAGER regards the processes of 
diffusion and flow of heat and electricity as being representable 
in terms of velocities, J;, which describe the rate at which 
the processes take place, and ‘thermodynamic forces,’ X;, 
which cause the processes to take place. He assumes that 
r 
linear relations, J; exist between the velocities 


and forces, the L;; being coupling coefficients. The theory 
bears a formal analogy to that of currents in electrical net- 
works, and the principle of microscopic reversibility yields 
jj=L;;. The X; 


the symmetrical relation L must be chosen 


ss 
so that . » J, X; is equal to the rate of generation of entropy 
| 

per unit volume. Chapter IV considers the rate of increase of 
entropy in various natural processes, and Chapter V deals 
with the application to a number of particular problems, such 
as thermal migration through a barrier, and thermal diffusion 
in electrolytic solutions. The final chapter discusses the 
relationship between the THomson hypothesis and the On- 
SAGER theory. Adequate references are given to original papers 
for those who wish to pursue the subject more thoroughly. 

The 


chemists, engineers and all others who are concerned with the 


book can be heartily recommended to physicists, 


applications of thermodynamics. C. Doms 


Combustion, 
Inc., 


Bernarp Lewis and GUENTHER ELBE; 
Flames and Explosions of Gases. 


Publishers, New York 1951. 


VON 
Academic Press 


795 pp. Price 5 14. 


This book is intended as a standard work on the physical 
and chemical principles of the combustion process. Its use is 
recommended to the student, engineer and research worker. 
It may be said that the authors have fully attained their aim. 
When comparing this work with the book of the same title 
published in 1938 and written by the same authors, it cannot 
properly be called a second edition, as the contents are entirely 
new, with the exception of a few minor sections, and the text 
has been considerably enlarged. This enlargement provides 
a good idea of the great activity in gas-combustion research 
and technique during the last ten years. 

Though a great many of the problems discussed in this 
book have not been solved definitely by far, the authors by 
no means give a dry summing-up of facts. Naturally, a large 
portion of the text is devoted to the experimental and theoretic- 
al work of the authors. However, due attention is also given 
to important work by others. Those who only want orientation 
in the field concerned will find in this book a well-founded 
representation of the most varied phenomena, problems and 
experimental methods which may be encountered in gas 
combustion. Those who want to enter more deeply into 
certain subjects will moreover find an extensive and up-to-date 
(1951) literature list. 
of numerical data is given in the form of tables and graphs. 


For practical workers a fine collection 


The arrangement of the text is largely the same as in the 
previous edition. The first part (about 200 pp.) is devoted to 
the reaction kinetics of the non-catalytic oxidation reactions 
H,, CO, hydrocarbons and some oxygen-containing organic 
substances are treated successively. In this part considerable 
attention is also paid to the important phenomenon called 
“engine knock.”’ 

The second part of the book deals with flame propagation 
(over 400 pp.). Much attention is given to the theory of the 
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flame front and to ignition. Moreover this part contains dis- 
cussions of diffusion flames, detonations and turbulence. 

The treatment of burners with air entrainment is of great 
practical importance. 

Part III deals with the state of the burned gas (60 pp.), 
while Part IV (30 pp.) is devoted to the technical application 
of gas combustion in industry and in engines, This last part 
is very brief and contains a survey of the basic principles only. 

The scope of the work under discussion does not allow of 
giving an extensive treatment of these subjects; this would 
require a separate book. 

The final 30 pp. are reserved for several valuable tables, 
containing ¢.g. thermodynamic data on gas reactions, flame 
temperatures and explosion limits. 

The production of the book is excellent. In spite of its 
price, which according to European standards is perhaps rather 
high, it should be bought by everyone concerned with the 


combustion of gases. Dr. J. VAN STEENIS 


Ullmanns Encyelopidie der Technischen Chemie, Dritte vollig 
neugestaltete Auflage (Editor: Dr. Witne tm Forrst) Band I 

Chemischer Apparatebau und Verfahrenstechnik, (Co-Editors: 
Prof. Dr. E. Wicker and Dipl.-Ing. E. ROémer). Publishers: 
Urban & Schwarzenberg, Miinchen-Berlin 1951. XXII 1011 
DM 10s 


pages, 18.5 « 26,5 em. 


Since the publi ation of the first edition of ULLMANNs 
famous handbook it has proved to be one of the great standard 
works in chemical literature. It is remarkable how even this 
first attempt was a success in every way; and since then this 
encyclopaedia has become indispensable for the industrial 
chemist and the chemical engineer as well as for the metall- 
urgical and mining engineer and workers in other fields of 
applied science. The last edition was published between 
1928 and 1932, so that at the moment the handbook is more 
than 20 vears old. It was reprinted after 1933 without revision, 
the only difference being that the name of the creator of 
this work was not allowed to appear on the title page. The 
work has found a successor in the U.S.: the Encyclopaedia of 


Kirk 


published from 1947 on with 13 volumes appearing at 7-month 


Chemical Technology, edited by and OreMerR and 
intervals. 


Owing to the impetuous development of chemical engineer- 
ing since 1930, ULLMANNS old Enevelopaedia is now out oft 


date in many respects, and for this reason the appearance of 
a new edition is very fortuitous. 

The first volume of the new edition has now been issued. 
Whereas the former editions of this handbook were arranged 
alphabetically according to catchwords, the first two of the 
new volumes show a different order. Volume I contains a sys- 
tematic survey of Chemical Engineering in a logical arrange- 
ment, Suecessive chapters deal with thermodynamics and 
flow of fluids, pressure and vacuum techniques, heating and 
cooling, separating and mixing processes; and finally a dis- 
cussion of reactors and another about construction materials 
are given. Naturally, many experts have contributed to the 
chapters devoted to their particular fields of study. 

Each chapter is preceded by a comprehensive table of 
contents, printed on thick vellow paper. Besides, two indexes 
have been inserted at the end of the book, one general and 
another for the chapter on construction materials. A list of 
authors is lacking, which is to be regretted, In many places 
literature references are given, naturally to a limited extent. 
European literature is given prevalent attention. 

Volume I is a beautiful book in a splendid binding, printed 
on high-quality paper and containing 1400 exceptionally fine 
figures. 

In all respects the first impression the book makes is a fav- 
ourable one. It is understandable that a work like this has 
its weeknesses. The treatment of thermodynamics is very 
brief, as is e.g. that of the separation of homogeneous gas 
mixtures (/.@. gas absorption). Presumably every specialist 
will hit on shortcomings in his own field, After all, with books 
of this scope the selection of the material to be treated is 
a matter of taste. Yet every one in possession of this work 
will find quick orientation possible. 

The book contains a number of splendid systematic and 
tabular surveys, ¢.g. of vacuum pumps, dimensionless heat 
transfer equations, crushing and grinding plants, etc. ChapterVI 
gives a complete and comprehensive survey of reaction appa- 
ratus, such as is not often met with elsewhere. 

To sum up, it should be said that editors and publishers 
may well be congratulated on the appearance of this book, 
which may be called indispensable for libraries. It is hoped 
that the volumes will appear at fairly short intervals and that 
they will live up to the standard promised by the first. How- 
ever, the price of works like this is so high that only very few 
private persons will be able to afford the luxury of buying 


them. LD. W. van KREVELEN 


ERRATA 
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IV by M.M. Harunt 


and J. ANDERSON STORROW. 


(1) in last six lines (column 2) page 157 
and first six lines (column 1) page 158 

“B” for “A” 

|] o's 

“A” for “C”. 


“ ubstitule 


(2) equation 9, page 161: 


include “4” in numerator of right hand side: 


to read “”qg” instead of “qq”. 





